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FOREWORD 


This  report  summarizes  the  work  accomplished  by  the 
Materials  Sciences  Corporation  under  NASA  Contract  NAS1-13931. 
Mr.  G.  L.  Roderick  (US  AAMRDL ,  Langley  Directorate)  was  the 
NASA  Project  Engineer. 

The  Program  Manager  and  Principal  Investigator  for 
Materials  Sciences  Corporation  were  Dr.  S.  V.  Kulkarni  and 
Dr.  R.  L.  Ramkumar,  respectively.  The  experimental  program 
was  conducted  at  the  University  of  Delaware  under  the  super¬ 
vision  of  Professor  R.  Byron  Pipes.  The  authors  would  like 
to  acknowledge  the  many  helpful  discussions  they  had  with 
Dr.  B.  W.  Rosen  and  Mr.  G.  L.  Roderick  during  this  program. 
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a 
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c 
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G  ■ 
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axial  shear  modulus  of  the  i  lamina  in 
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unnotched  laminate  tensile  strength 

unnotched  laminate  compressive  strength 
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T 

xyi 


yield  value  of  shear  stress  in  the  i 
lamina  in  the  xy  coordinate  system 


txz12 ' Txz23 


yield  value  of  interlaminar  shear  stress 
between  layers  1  and  2,  and  2  and  3, 
respectively. 


EVALUATION  AND  EXPANSION  OF  AN  ANALYTICAL  MODEL 
FOR  FATIGUE  OF  NOTCHED  COMPOSITE  LAMINATES 

R.  L.  Ramkumar,  S.  V.  Kulkarni* 
and  R.  B.  Pipes** 

Materials  Sciences  Corporation 
SUMMARY 

This  report  describes  the  analytical  and  experimental 
study  performed  to  expand  the  existing  static  and  fatigue 
failure  analysis.  The  analytical  effort  extended  the  analysis 
to  include  interlaminar  effects,  while  the  experimental  effort 
developed  methods  to  obtain  basic  experimental  data  required 
as  input  to  the  analysis. 

The  static  failure  analysis  for  notched  laminates  was 
modified  to  include  interlaminar  effects  near  the  notch.  The 
interlaminar  effects  were  included  by  discretizing  the  laminate 
into  laminae  in  the  region  of  the  projection  of  the  notch  width 
along  the  axial  direction.  This  discretization  produced  addi¬ 
tional  degrees  of  freedom  which  were  used  to  analyze  axial 
in-plane  shear  damage  or  crack  in  each  lamina  and  axial  inter¬ 
laminar  damage  or  delamination  at  the  laminae  interfaces. 

Because  a  large  number  of  damage  zones  is  possible  when  elastic- 
perfectly  plastic  behavior  in  in-plane  and  interlaminar  shear 
is  considered,  the  analysis  was  simplified  by  assuming  a 
secant  modulus  representation  of  the  nonlinear  shear  stress- 
shear  strain  curve.  The  growth  of  damage  by  both  the  above- 
mentioned  approaches  is  calculated  for  various  realistic 
T-300  Gr/Ep  laminates  (  [0~/±45]  ,  [0/±45/90]  )  and  compared 

with  earlier  experimental  observations. 

Three-dimensional  elastic  and  two-dimensional  elastic- 
plastic  finite  element  analyses  were  also  performed  for  some 
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**Associate  Professor,  University  of  Delaware,  and  Affiliate 
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notched  laminates  for  the  purposes  of  (i)  developing  an  under¬ 
standing  of  the  behavior  of  interlaminar  stresses  around  the 
notch,  and  (ii)  determining  the  nature  of  the  2-D  elastic- 
plastic  stress  distribution  adjacent  to  the  notch  when  the 
maximum  attainable  stresses  or  combinations  thereof  for  the 


laminate  are  governed  by  a  yield  criterion.  The  results 
obtained  from  task  (i)  provided  a  guideline  for  qualitatively 
assessing  the  approximate  stress  analysis  procedure.  Task 

(ii)  determined  how  appropriate  some  of  the  postulates  (such 
as  the  average  stress  concentration  region  being  a  function  of 
of  laminate  anisotropy)  of  the  basic  failure  model  are. 

The  experimental  effort  was  directed  at  several  differ¬ 
ent  tasks:  (i)  determining  static/fatigue  (S-N)  curve  data 

for  interlaminar  shear  and  normal  stresses;  (ii)  assessing 
the  effect  of  (a)  the  interaction  between  in-plane  shear  and 
transverse  compression  on  the  wearout  rates  for  strength  and 
stiffness  in  a  [ ± 4 5 ) s  laminate  subjected  to  uniaxial  loading, 
and  (b)  the  transverse  restraint  (as  provided  by  90°  layers 
in  a  [0i/90.]g  laminate)  on  the  failure  mode  and  wearout  rates 
for  strength  and  stiffness  for  a  unidirectional  layer  when 
subjected  to  compression/compression  fatigue;  and  finally 

(iii)  designing  the  experiments  and  specimens  to  obtain  the 
required  data. 
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INTRODUCTION 


The  use  of  composite  materials  for  primary  aerospace 
structural  applications  requires  that  they  have  adequate  safety 
margins  and  lifetimes  under  all  anticipated  environments. 
Accordingly,  a  sound  understanding  of  the  fracture  and  fatigue 
behavior  of  these  materials  is  needed,  and  the  materials  must 
also  be  well  characterized  by  soundly  based  experimental  data. 
To  characterize  materials  by  tests  alone  is  expensive;  analy¬ 
ses  which  correlate  successfully  among  data  points  can  greatly 
reduce  the  number  of  tests  needed  and  guide  the  testing 
methods  to  provide  more  useful  experimental  results. 

The  development  of  such  analyses  was  initiated  in  refer¬ 
ence  1  in  which  a  methodology  was  developed  for  predicting 
the  growth  of  cracks  and  ultimate  failure  under  fatigue  loading 
of  fiber  composite  laminates  containing  through-the- thickness 
notches  or  holes.  This  framework  was  based  on  knowledge  of 
the  static  and  fatigue  behavior  of  a  unidirectional  layer  of 
the  material.  Specifically,  the  analysis  assumes  that  lamina 
fatigue  behavior  within  the  notched  laminate  is  predictable 
from  fatigue  data  obtained  from  tests  of  unnotched  laminae. 

To  verify  this  analysis,  and  develop  a  semi-empirical 
capability  for  the  fatigue  analysis  of  notched,  fiber-composite 
laminates,  the  predictions  derived  from  reference  1  were  com¬ 
pared  in  reference  2  with  experimental  fatigue  data  for  notched 
[02/±45]g  Boron/epoxy  laminates  containing  a  0.635  cm  diameter 
circular  hole.  The  specific  phenomena  studied  were: 

(1)  initiation  of  fatigue  damage,  and  growth  as  a  func¬ 
tion  of  load  cycles  ; 

(2)  fatigue  life  and  mode  of  failure;  and 

(3)  residual  strength  and  mode  of  failure  after  a  pre¬ 
determined  number  of  cycles. 

The  residual  strength  after  fatigue  loading,  and  both 
axial  and  transverse  damage  growth,  were  measured  quantita¬ 
tively  for  correlation  with  the  analytical  predictions. 
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Correlation  between  analysis  and  experiment  was  hampered, 
however,  by  the  following  factors: 

(1)  delamination  of  the  0°  surface  layers  in  the  region 
of  the  longitudinal  projection  of  the  notch  diameter 

(2)  lack  of  a  statistically  significant  data  base  for 
lamina  fatigue  properties  in  general;  and,  par¬ 
ticularly  , 

(3)  absence  of  lamina  axial  compression  fatigue  data. 

The  principal  deficiency  in  the  analysis  appeared  to  be 
the  lack  of  a  capability  to  predict  the  growth  of  the  delamina 
tions.  A  primary  objective  of  the  present  study  is  to  correct 
this  deficiency.  The  approach  used  is  to  modify  the  existing 
static  failure  analysis  (ref.  3)  to  permit  the  effects  of 
interlaminar  shear  stress  to  be  taken  into  account.  This  is 
achieved  by  discretizing  the  laminate  into  laminae  in  the 
vicinity  of  the  notch.  Significant  features  of  the  analytical 
developments  are  described  in  the  section  '‘Description  of  the 
Failure  Model,"  and  the  mathematical  details  are  presented  in 
Appendix  A.  The  analysis  has  been  incorporated  into  a  com¬ 
puter  program  for  calculation  of  damage  growth  and  residual 
strength. 

To  assess  the  degree  of  realism  associated  with  the 
assumptions  in  the  analysis,  the  calculated  stresses  have 
been  compared  to  those  obtained  from  2-D  and  3-D  finite  ele¬ 
ment  analyses.  These  comparisons  are  discussed  in  the  sec¬ 
tion  "Finite  Element  Analysis  of  Notched  Composite  Laminates." 

An  experimental  program  was  conducted  to  demonstrate 
the  feasibility  of  obtaining  the  data  required  as  input  to 
the  analysis.  Inadequacies  of  the  existing  data  base  for 
lamina  fatigue  properties  (already  noted  during  the  analysis- 
experiment  correlation  study  in  reference  2)  are  compounded, 
by  the  need  to  incorporate  consideration  of  interlaminar 
effects  into  the  fatigue  analysis.  Exploratory  testing  to 
develop  methods  for  obtaining  lamina  properties  from  simple 


4 


laminates  instead  of  unidirectional  samples  therefore  has  been 
initiated  as  follows: 


(i)  [±45]  laminate  uniaxial  compression/compression 

fatigue  tests  have  been  made  to  determine  the  nature 
of  the  interaction  between  in-plane  shear  and  trans¬ 
verse  compression  in  a  lamina  and  to  compare  the 
wearout  rate  with  that  of  a  [±45]g  uniaxial  tension/ 
tension  fatigue  test; 


(ii)  [CL/90j]g  laminate  uniaxial  compression/compression 
fatigue  tests  have  been  made  to  ascertain  the  effect 
of  transverse  restraint  (provided  by  the  90-degree 
laminae)  on  the  failure  mode  and  the  wearout  rate 
for  a  unidirectional  lamina  in  compression/ 
compression  fatigue; 


(iii)  [+252/90]  laminate  specimens  were  used  to  obtain 

the  S-N  curve  for  interlaminar  normal  tensile  stress; 
and 


(iv)  Partially  notched  [±45./0-]  laminates  were  used 

1  D  s 

for  determining  the  behavior  of  interlaminar  shear 
stress  in  fatigue. 

Highlights  of  the  tests  and  their  implications  in  the 
light  of  the  analysis  as  extended  herein  are  discussed  in  the 
section  entitled  "Experimental  Program."  The  test  specimen 
design  is  outlined  in  Appendix  B. 
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FAILURE  MODEL  FOR  NOTCHED  COMPOSITE  LAMINATES 
INCLUDING  INTERLAMINAR  EFFECTS 


DISCUSSION  OF  PREVIOUS  MODEL 

The  prior  analysis  (ref.  2,  see  fig.  1)  was  based  on  two 
concepts:  (1)  Failure  in  the  presence  of  a  notch  can  result 

from  cracks  propagating  along  planes  of  weakness  parallel  to 
the  fibers  or  perpendicular  to  them;  (2)  The  notched  laminate 
may  be  divided  into  regions  having  simple,  approximate  stress 
states.  The  first  of  these  concepts  is  utilized  unchanged  in 
this  study.  The  second  is  extended  to  take  into  account 
through-the-thickness  effects. 

In  reference  2,  the  laminate  was  considered  as  an  assem¬ 
blage  of  discrete  regions;  namely,  the  strip  containing  the 
notch;  the  overstress  region  adjacent  to  this  strip;  the  inter¬ 
mediate,  shear-overstress  region;  and  the  undisturbed  or  aver¬ 
age  stress  region.  In  the  thickness  direction,  the  laminate 
was  treated  as  a  homogeneous  material.  Consequently,  the  pre¬ 
dictions  of  interlaminar  stresses  and  the  resulting  delamina¬ 
tions  were  beyond  the  scope  of  the  analysis.  To  take  into 
account  the  interlaminar  stresses,  the  laminate  was  considered 
herein  divided  into  laminae  in  the  local  regions  where  the 
effect  of  interlaminar  stresses  is  most  pronounced.  Resultant 
changes  in  the  analysis  and  its  failure-prediction  capability 
are  discussed  in  the  following  section. 

DESCRIPTION  OF  THE  MEW  FAILURE  MODEL 

Static  Fa. ilure  Analysis 

The  static  failure  model,  at  the  laminate  level,  is 
divided,  as  before  (ref.  1) ,  into  three  regions:  (i)  a  central 

core  region  which  is  the  projection  of  the  notch  in  the  loading 
direction;  (ii)  an  overstressed  region  of  average  stress  con¬ 
centration,  adjacent  to  the  core  region;  and  (iii)  an  average 
stress  region  (fig.  1) .  As  mentioned  in  the  earlier  sections, 
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than  two  damage  zones  exist  in  the  laminate.  However,  this 
is  alleviated  by  eliminating  a  few  damage  sequences  that  are 
unlikely  to  occur.  It  is  assumed  that,  before  too  many  damage 
zones  appear,  the  load  will  be  large  enough  to  cause  a  trans¬ 
verse  failure  or  one  of  the  earlier  defects  will  precipitate 
an  axial  failure.  This  is  where  the  simplicity  in  the  applica¬ 
tion  of  the  elastic,  secant  modulus  approach  is  obvious.  The 
elimination  of  the  inelastic  regions  leaves  only  cracks  and 
delaminations  to  be  accounted  for. 

In  many  laminate  constructions ,  a  degenerate  region  is 
created  whenever  the  applied  loading  reaches  a  critical  value. 
Degeneracy  sets  in  when  one  or  more  layers  in  the  core  region 
are  surrounded  by  stress-free  surfaces  and/or  inelastic  sur¬ 
faces.  In  the  elastic-plastic  approach,  one  or  more  degenerate 
layers  have  a  linear  variation  of  axial  strain  in  the  degenerate 
region.  The  degenerate  layers  in  the  elastic,  secant  modulus 
approach  are  stress  free  and  are  essentially  "carried"  by 
adjacent  load-carrying  regions.  In  a  degenerate  region,  the 
solutions  for  the  degenerate  layers  are  uncoupled  from  the 
five  basic  governing  equations.  The  solutions  change  for  dif¬ 
ferent  types  of  degeneracies  (Appendix  A) .  Hence,  if  all  types 
of  laminate  behavior  are  to  be  included  in  the  solution  pro¬ 
cedure,  programming  of  the  solution  becomes  tedious. 

Fatigue  Failure  Analysis 

A  "mechanistic  wearout"  concept  used  in  earlier  studies 
(ref.  2)  underlines  the  basic  philosophy  of  the  laminate 
fatigue  behavior.  On  fatigue  loading,  material  property  degra¬ 
dation  is  predominant  in  the  vicinity  of  the  notch  due  to 
stress  concentration  effects.  An  experimental  characteriza¬ 
tion  of  the  lamina  fatigue  data  helps  estimate  the  degrada¬ 
tion  in  the  laminate  properties  due  to  cyclic  loading.  The 
material  property  degradation  on  fatigue  loading,  when  incor¬ 
porated  into  the  static  failure  model  described  earlier,  could 
lead  to  fatigue  failure  modes  and  strengths  that  are  very 
different  from  the  static  predictions. 
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Figure  6  shows  a  flow  chart  of  the  program  developed  to 
predict  the  failure  modes  and  strengths  of  notched  laminates 
subjected  to  cyclic  loading. 

DISCUSSION  OF  RESULTS 

The  approximate  failure  analysis  for  notched  laminates 
which  considers  the  effects  of  interlaminar  stresses  was 
utilized  to  investigate  the  damage  growth  in  some  T-300/5208 
laminates.  Table  1  lists  the  various  data  sets  utilized  in 
the  failure  analysis.  The  trial  data  are  used  to  illustrate 
the  extent  of  the  predictive  capability  of  the  analysis  when 
laminate  failure  in  the  transverse  direction  is  prevented. 

Tables  2  and  3  indicate  the  results  of  the  analysis  for 
T-300/5208  [02/±45/0]s  and  [±45/02/0]  laminates,  respectively, 

for  different  values  of  interlaminar  boundary  layer  zone  size 
"d"  and  the  trial  data.  In  the  table,  the  distance  "d"  is 
the  interlaminar  boundary  layer  region  which  assumes  two 
values:  (i)  half  the  lamina  thickness  (typical  for  Boron/ 

epoxy);  and  (ii)  one-tenth  the  lamina  thickness  (typical  for 
Graphite/epoxy) .  Figure  7  shows  a  typical  variation  of  the 
in-plane  displacement  as  obtained  in  reference  4.  The  figure 
explains  the  meaning  of  the  parameter  "d".  In  the  present 
analysis,  the  axial  displacements  in  the  various  layers  are 
piecewise  uniform,  and  hence,  discontinuous  at  the  interfaces. 
The  ratio  of  the  displacement  discontinuity  to  "d"  -  the  trans¬ 
verse  shear  strain  -  should  match  the  slope  of  the  continuous 
displacement  curve  at  the  interface  obtained  from  a  higher 
order  analysis.  For  the  [C^/ilS/O^  laminate,  the  in-plane 
plastic  zone  is  initiated  almost  simultaneously  in  the  C>2  and 
0  laminae  (a  ^>0)  •  Also,  the  next  damage  zone  -  a  plastic 
zone  -  appears  in  the  ±45  laminae  (a  2^)  because  the  inter¬ 
laminar  shear  stress  (t^  n„)  is  not  large  enough  to  form  an 
interlaminar  inelastic  region.  Thus,  the  approximation  that 
the  axial,  in-plane  plastic  zone  grows  uniformly  through- the- 
thickness  (ref.  3)  may  reflect  the  actual  behavior  for  this 
laminate.  After  both  the  plastic  zones  have  grown  over  a 
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certain  distance,  an  in-plane  crack  is  initiated  in  the  C>2 
and  0  layers.  Note,  however,  that  the  interlaminar  shear 
stresses  are  high  and  may  eventually  give  rise  to  a  delamina¬ 
tion.  In  tables  2  and  3,  the  results  for  the  [02/±45/0]s 
laminate  are  essentially  the  same  as  for  the  [±45/02/0]s 
laminate  except  that  the  magnitude  of  the  interlaminar  shear 
stress  (txz]_2^  s  increased  by  about  25%. 

An  important  conclusion  which  can  be  drawn  from  the  above 
results  is  that  the  number  of  combinations  of  different  damage 
zones  (inelastic  regions  as  well  as  cracks)  is  large  (fig.  5) . 
The  factors  which  affect  this  number  are  the  material  proper¬ 
ties,  stacking  sequence,  and  notch  size.  An  automated  numer¬ 
ical  scheme  to  "track"  the  various  possible  damage  growth 
combinations  is  a  tedious  task.  Also,  the  solution  procedure 
expands  with  the  initiation  of  plastic  zones.  Thus  the  analy¬ 
sis  is  further  complicated.  Consequently,  the  analytical 
procedure  was  modified  to  eliminate  the  need  to  consider 
plastic  zones  by  considering  linear  approximations  of  the 
nonlinear  in-plane  and  inter laminar  shear  stress-strain  curves 
as  shown  in  figure  8.  Such  an  analysis  would,  however,  still 
consider  the  growth  of  in-plane  shear  cracks  and  interlaminar 
debonds,  but  the  number  of  combinations  of  the  different 
damage  zones  is  drastically  reduced.  The  various  damage  growth 
combinations  are  illustrated  in  figure  5. 

Figures  9  and  10  illustrate  the  progressive  damage  growth 

for  [0_/±45/0]  and  [0/±45/90]  T-300/5208  Graphite/epoxy 

Z  s  s 

laminates,  respectively,  based  on  the  trial  data  in  table  1. 

The  force  represents  the  applied  laminate  force  and  "a" 
indicates  the  cumulative  sum  of  the  damage  zones.  In  figure 
9,  two  inelastic  in-plane  axial  damage  zones  are  predicted 
before  an  axial  crack  in  the  0°  layers  appears.  Figure  10 
illustrates  a  similar  pattern  of  damage  growth,  except  that 
no  in-plane  crack  occurs  in  the  90°  layer  and  an  interlaminar 
inelastic  zone  is  initiated  at  a  higher  value  of  the  applied 
force.  Figure  11  shows  the  effect  of  notch  size  on  the  damage 
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growth  in  a  [±45/0/0]  laminate  for  the  trial  data  input.  The 

[±45/0/0]  notation  is  used  for  a  [±45/0-]  laminate,  which 

z  s 

is  treated  by  the  analysis  as  a  "±45°  layer"  and  two  0°  layers 
If  the  ±45  layup  is  treated  as  two  layers  (a  +45  layer  and 
a  -45  layer)  in  the  analysis,  the  Poisson  effect  causing  an 
increase  in  the  shear  modulus  will  not  be  accounted  for. 
Figures  9  to  11  show  the  capability  of  the  approximate  analysi 
procedure  to  predict  damage  regions  on  a  lamina  level.  How¬ 
ever,  as  has  been  noted  earlier,  the  analysis  is  not  capable 
of  predicting  all  combinations  of  damage  growth  because  of  a 
significant  increase  in  the  complexity  of  the  computation 
procedure . 

The  variation  of  failure  load  with  notch  size  is  studied 

in  figure  12  for  a  [±45/0/0]  laminate  for  the  actual  data 

s 

input.  The  following  conclusions  can  be  drawn  from  the  figure 

(i)  Damage  zone  size  increases  with  increase  in  notch 
size  . 

(ii)  In-plane  damage  occurs  only  in  the  0°  layers. 

(iii)  Laminate  fails  by  transverse  crack  propagation. 

(iv)  At  failure,  the  maximum  in-plane  shear  stress  in 
the  ±45  laminae  is  higher  for  larger  notch  sizes 
than  that  for  smaller  notch  sizes,  thus  indicating 
the  increased  probability  of  through-the-thickness 
axial  crack  propagation  for  larger  notch  sizes. 

(v)  The  transverse  failure  values  predicted  by  the 

secant  modulus  approach  are  close  to  those  predic¬ 
ted  by  the  elastic-plastic  analysis. 

Figure  13  illustrates  the  results  for  a  [0/0/ ±45] 

s 

laminate  with  the  allowable  interlaminar  shear  stress  value 
reduced  to  account  for  the  presence  of  interlaminar  normal 
tension  stress  at  the  0/±45  interface.  The  obvious  result 
is  the  initiation  of  an  interlaminar  inelastic  zone  followed 
by  in-plane  plasticity  in  the  0°  layers. 
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In  figure  14  the  results  for  a  [0/±45/90]g  laminate  for 
different  notch  sizes  are  shown  for  the  actual  data.  Since 
this  laminate  is  more  isotropic  than  the  [ 0 2/ ±4 5 ] s  laminate, 
the  inelastic  damage  zones  are  smaller.  Also,  failure  occurs 
by  transverse  crack  propagation. 

Figure  14  is  an  illustration  of  the  predictive  capability 
of  the  secant  modulus  analysis.  As  shown  in  the  previous 
figure,  the  [0/±45/90]s  laminate  fails  in  the  transverse 
direction  before  any  crack  appears  when  the  secant  modulus 
approach  was  used.  In  figure  15,  the  damages  are  made  to 
grow  by  preventing  a  transverse  failure.  Table  4  shows  the 
effect  of  the  interlaminar  boundary  layer  thickness  on  the 
damage  growth  for  the  same  laminate. 

Figure  16  illustrates  an  interesting  interlaminar  damage 
growth  in  a  [0/0/±45]  laminate  for  a  reduced  value  of  inter¬ 
laminar  stress  that  causes  delamination.  The  reduction 
accounts  for  the  effect  of  interlaminar  normal  tensile  stress. 
A  high  value  of  a  ,  the  applied  stress  at  which  transverse 
failure  occurs,  is  chosen  to  prevent  that  failure  mode.  It 
is  seen  that  the  delamination  that  initiates  at  the  0/±45 
interface  grows  in  an  unstable  fashion  to  failure.  The  0° 
plies  thus  separate  away  from  the  ±45  plies  in  the  notch 
region.  This  result  was  observed  experimentally  (fig.  17) 
in  the  first  phase  of  the  program  (ref.  2) .  This  qualitative 
agreement  between  the  failure  prediction  of  the  analytical 
model  and  the  experimental  observation  justifies  the  approach 
taken  in  building  the  model. 

In  the  results  discussed  above,  the  ultimate  failure 
stress  for  the  ±45  laminate  has  been  obtained  from  "netting 
analysis."  The  netting  analysis  predictions  may  be  sig¬ 
nificantly  higher  than  first-ply  failure  predictions  because 
fiber  failure  within  a  ply  is  required  in  "netting  analysis." 
Also,  a  through-the- thickness  axial  crack  cannot  propagate 
in  a  ±45  laminate  unless  the  fibers  fail.  Laminate  failure 
stress  corresponding  to  the  fracture  of  ±45  layer  fibers  is 
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synonymous  to  the  netting  analysis  failure  stress.  Hence,  the 
upper-bound  predictions  of  the  netting  analysis  may  be  realis¬ 
tic  when  considering  the  incipient  conditions  required  for 
the  initiation  and  growth  of  through- the-thickness  axial  cracks. 
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FINITE  ELEMENT  ANALYSIS  OF 
NOTCHED  COMPOSITE  LAMINATES 


The  failure  predictive  capability  of  the  analytical  model 
was  based  on  a  few  simplifying  assumptions.  Justification  of 
these  assumptions  is  essential  to  extend  the  model  application 
to  any  general  laminate. 

A  three-dimensional  finite  element  analysis  was  carried 
out  to  predict  the  in-plane  and  interlaminar  stresses  in  a 
notched  laminate.  The  primary  objectives  of  the  3-D  finite 
element  analysis  were  (i)  to  substantiate  the  assumptions 
made  in  the  approximate  model  and  (ii)  to  provide  guidance 
for  any  modification.  This  elastic  analysis  aimed  at  estab¬ 
lishing  a  qualitative  agreement  with  the  model  stress  pre¬ 
dictions,  and  also  establishing  the  predominance  of  the  axial 
normal  and  interlaminar  shear  stresses  over  the  other  stresses. 
A  two-dimensional,  elastic-perf ectly  plastic  finite  element 
analysis  was  carried  out  to  verify  the  inelastic  damage  growth 
predictions  of  the  analytical  model  for  various  notched 
laminates . 

THREE-DIMENSIONAL  ELASTIC  ANALYSIS 

The  present  analytical  modeling  philosophy  considers  the 
physics  of  the  problem  without  accounting  for  the  detailed 
stress  distribution.  Thus,  only  those  stress  fields  (and  the 
appropriate  equations  of  equilibrium)  are  considered  which 
have  a  predominant  effect  on  the  observed  failure  modes.  The 
driving  force  for  this  approximate  approach  has  been  the  com¬ 
plexity  of  the  three-dimensional  nonlinear  stress  analysis  and 
the  associated  failure  prediction.  However,  the  use  of  simpli¬ 
fied  models  (which  neglect  the  effects  of  some  stress  fields) 
can  be  justified  only  if  there  exists  qualitative,  and  to  a 
limited  extent,  quantitative  correlation  with  the  more  exact 
stress  analyses.  To  that  end,  a  three-dimensional  finite 
element  anlaysis  of  notched  laminates  was  carried  out  to  check 
the  assumptions  made  in  the  approximate  stress  analysis  pro¬ 
cedure. 
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The  computer  code  SAP-IV  was  used  to  carry  out  the  finite 
element  analysis  (ref.  5) .  An  octant  of  the  notched  laminate 
was  considered  and  each  layer  above  the  reference  midplane 
was  discretized  into  the  mesh  form  shown  in  figure  18 .  A 
minimum  of  eight  and  a  maximum  of  twenty-one  nodes  were  chosen 
to  describe  the  general  three-dimensional  isoparametric  ele¬ 
ment  used  in  the  analysis.  The  element  was  used  to  represent 
orthotropic,  elastic  media.  Laminate  analysis  was  used  to 
apply  the  appropriate  loads  in  each  layer  so  that  a  far-field 
uniform  strain  state  existed  in  the  laminate. 


In  figure  19,  a  obtained  from  the  approximate  analysis, 
y  z 

at  the  0/±45  interface,  is  compared  with  the  a  variation 

y  z 

from  the  3-D  finite  element  analysis  (the  average  value  at 
the  face  centroids) .  The  comparison  is  made  at  two  axial 
locations.  As  shown  in  figure  19,  the  average  stresses  from 
the  approximate  analysis  compare  qualitatively  with  the  3-D 
finite  element  results. 

Figure  20  shows  the  axial  normal  stress  variation  in  the 
load-carrying  direction  agreeing  qualitatively  with  the  analyti¬ 
cal  prediction.  Figure  21  shows  the  variation  of  the  same 
stress,  a  ,  in  the  elements  just  above  the  x  axis.  A  com¬ 
parison  with  the  applied  stress  in  each  layer  (fig.  18)  gives 
a  measure  of  the  stress  concentration  in  each  of  the  three 
layers  of  the  [0/+45/-45]s  laminate.  These  stress  concentra¬ 
tion  values  show  good  agreement  with  the  value  predicted  by 
the  analytical  model. 

The  number  of  integration  points  to  be  used  in  the  numeri¬ 
cal  evaluation  of  integrals  over  volumes  is  referred  to  as  the 
integration  order  in  SAP-IV.  Different  orders  of  integration 
may  be  chosen  in  the  plane  of  the  element  and  through  its 
thickness.  The  orders  of  integration  chosen  in  the  finite 
element  analysis  affect  the  cost  of  the  run  and  the  magnitude 
of  the  results.  Figures  22  and  23  show  that  for  two  different 
orders  of  integration  (2  and  4) ,  the  0 yZ  variations  with  y  at 
the  0/45  interface  for  a  [0/±45/-45]  laminate  have  good  quali¬ 
tative  agreement  along  with  a  quantitative  difference. 
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The  approximate  analysis,  discussed  in  the  previous  chap¬ 
ter,  treats  notches  of  any  geometry  as  equivalent  slit  notches 
since  the  notched  strength  for  composite  laminates  is  not 
strongly  dependent  on  the  notch  shape  (ref.  3).  The  3-D 
finite  element  analysis  was  used  to  verify  this.  Figure  24 
compares  the  shear  stress  variations  for  a  circular  hole  and 
a  slit  notch  of  the  same  width  at  the  ±45/-45  interface  of 
the  [0/±45/-45]  laminate.  The  results  show  a  slight  differ¬ 
ence  only  in  the  vicinity  of  the  notch,  justifying  the  notch 
geometry  assumption  in  the  simplified  analytical  model. 

The  3-D  finite  element  results  indicate  that  the  variation 
of  in-plane  stresses  (a  an<3  aXy)  justifies  the  assumption 
of  an  average  stress  concentration  region  and  a  localized 
region  of  shear  stress  transfer.  The  finite  element  solutions 
also  justify  the  approximate  analysis  assumption  of  a  uniform 
value  of  °yZ  (independent  of  the  x-coordinate  but  a  function 
of  the  y-coordinate)  in  the  notched  region  and  a  zero  value 
in  the  other  regions  of  the  laminate.  The  non-zero  value  of 
a  predicted  by  the  3-D  finite  element  analysis  in  the  region 
away  from  the  notch  is  of  little  significance  because  it  does 
not  contribute  to  the  precipitation  of  failure. 

2-D  ELASTIC-PLASTIC  FINITE  ELEMENT  ANALYSIS 

As  an  additional  task  conducted  during  the  present  study, 
a  2-D  elastic-plastic  finite  element  analysis  was  performed 
for  some  notched  laminates.  Because  of  the  consideration  of 
elastic-perf ectly  plastic  behavior  in  shear  in  the  basic 
fracture  model,  the  elastic-plastic  finite  element  analysis 
was  used  to  investigate  the  damage  growth  and  stress  distribu¬ 
tion  near  the  notch  and  to  ascertain  the  effect  of  laminate 
anisotropy  on  damage  propagation.  The  PLANE  finite  element 
program  (developed  at  Grumman,  ref.  6)  was  used  for  this 
purpose . 

Figure  25  illustrates  the  finite  element  mesh  and  the 
loading  for  a  laminate  containing  a  slit  notch.  Input  values 
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of  moduli  and  strength  (ref.  7)  for  various  laminates  are  shown 
in  table  5.  Note  that  the  failure  stresses  satisfy  the  stabil¬ 
ity  condition  for  the  Hill's  yield  criterion  (ref.  8) .  Results 
of  the  finite  element  analysis  are  presented  for  a  [0/±45]g 
laminate  in  figures  26  and  27.  The  notch  blunting  effect  due 
to  plasticity  is  evident  in  both  the  figures.  The  figures  also 
show  the  change  in  the  notch  blunting  effect  on  varying  the 
applied  load.  Subsequently,  the  plastic  region  growth  is  inves¬ 
tigated  for  laminates  with  various  anisotropies.  The  critical 
value  of  the  loads  (a  multiple  for  the  applied  load  vectors)  at 
which  plasticity  is  initiated  for  the  [0],  [0/90] s  and  [02/±45]s 

T-300/5208  laminates  is  tabulated  in  table  6.  As  expected, 
the  load  is  lower  for  the  [0]  and  [0/90]  laminates.  In 
figures  28  through  30  are  shown  the  elements  which  have  yielded 
for  all  the  three  laminates  for  similar  load  levels.  It  is 
evident  that  the  yield  zones  are  most  prominent  for  the  [0/90] 
and  [0]  laminates,  in  that  order.  However,  it  has  been  ob¬ 
served  in  experiments  with  Boron/epoxy  laminates  that  the 
axial  inelastic  damage  zone  is  more  pronounced  in  the  [0] 
laminates  than  in  the  [0/90]  laminate  (ref.  9) .  This  is 
also  predicted  by  the  static  failure  model,  when  the  effec¬ 
tive  in-plane  shear  modulus  is  modified  for  laminates  contain¬ 
ing  90°  layers  in  order  to  reflect  the  local  bending  effects 
of  90°  fibers  in  the  vicinity  of  the  notch  (ref.  3) .  The  2-D 
finite  element  analysis  does  not  account  for  these  changes. 
Hence,  it  appears  that  the  elastic-plastic  finite  element 
analysis  is  not  capable  of  predicting  accurately  the  growth 
of  damage  zone  for  the  [0]  laminate.  Table  7  shows  the  varia¬ 
tion  of  an  approximate  measure  of  stress  concentration  factor 
with  increasing  load  for  different  laminates.  The  stress  con¬ 
centration  shows  a  decrease  with  increasing  load  for  the 
[0/90]  and  [02/±45]c  laminate,  while,  for  the  [0]  laminate, 
it  increases.  This  increase  is  contrary  to  the  well-known 
notch  blunting  effect  when  there  is  a  growth  of  the  inelastic 
region  in  the  axial  direction. 
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EXPERIMENTAL  PROGRAM 


The  experimental  tasks  conducted  during  this  program  have 

been  defined  by  the  observations  made  during  the  previous 

analysis-experiment  correlation  study  and  the  expansion  of  the 

basic  failure  model  to  include  interlaminar  effects.  The 

tests  are  of  an  exploratory  nature  and  were  performed  for  the 

specific  purpose  of  defining  the  requirements  of  basic  data 

to  be  utilized  as  input  in  the  analysis. 

The  thrust  of  the  earlier  experimental  program  (ref.  2) 

was  to  provide  Boron/epoxy  lamina  static/fatigue  data  to  be 

used  as  input  into  the  fatigue  analysis  and  to  generate  static/ 

fatigue  data  for  notched  Boron/epoxy  laminates  to  be  utilized 

in  the  fatigue  failure  model.  However,  in  the  previous  programs 

several  items  were  not  investigated. 

Consequently,  there  is  a  need  for  uniaxial  compression 

fatigue  data,  even  for  a  tension/tension  fatigue  loading  of  a 

notched  laminate.  To  that  end,  [0]  coupon  compression  tests 

were  performed.  However,  since  failure  in  a  [0]  longitudinal 

compression  test  is  accelerated  by  the  presence  of  longitudinal 

cracks,  the  residual  strength/stiffness  data  obtained  may  not 

be  representative  of  a  unidirectional  layer  in  a  multi-layered 

laminate,  because  of  the  absence  of  the  restraining  effect 

provided  by  the  adjacent  layers.  Hence,  utilization  of  the 

unidirectional  compression  data  alone  to  predict  laminate 

fatigue  response  is  questionable.  Alterations  in  failure  modes 

for  static  compression  have  been  observed  in  reference  10  on 

providing  a  nominal  lateral  reinforcement.  For  these  reasons, 

it  may  be  imperative  to  obtain  the  compressive  fatigue  data 

in  the  presence  of  transverse  restraint.  Indeed,  just  as  the 

[±45]  specimens  are  used  for  shear  data,  another  simple 
s 

laminate  may  be  required  for  other  layer  properties.  For  axial 

compression,  the  [0./90.]  laminate  provided  a  means  to  do  this. 

i  3  s 

Comparison  of  [0./90.]  and  [0]  laminate  compression  fatigue 

i  1  s 

test  results  will  indicate  the  effect  of  lateral  constraint  on 
residual  properties.  (Note  that  the  residual  properties  for  the 
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0°  layers  in  the  [0^/90 j]  specimen  have  to  be  extracted  by 
filtering  out  the  contribution  of  90°  layers) . 

Combined  stress  effects  upon  wearout  have  not  been  evalu¬ 
ated  experimentally  and  have  been  modeled  approximately  in 
the  fatigue  model.  Hence,  initial  assessment  of  the  magnitude 
of  this  problem  was  obtained  from  compression  fatigue  tests 
on  [±45]  laminates.  These  results  can  be  compared  with  those 
obtained  from  tension  fatigue  tests  on  similar  laminates.  Thus 
axial  in-plane  shear  combined  with  transverse  tensile  stress 
will  be  compared  with  the  same  shear  and  transverse  compressive 
stress . 

Important  parameters  in  the  analysis  which  considers  debond 
propagation  are  the  interlaminar  strengths  (shear  and  normal 
tension) .  An  estimate  of  the  interlaminar  shear  strength  may 
be  obtained  from  the  short  beam  test.  Alternatively,  the 
interlaminar  shear  strength  was  obtained  from  a  [±45/0]  lam¬ 
inate  with  a  slit  (transverse  to  the  fibers)  in  the  0°  layers. 
Also,  as  suggested  in  reference  11,  the  interlaminar  (a  ) 
normal  tension  strength  was  determined  by  utilizing  a  [±252/90] 
laminate.  The  [±252/90]  laminate  is  optimized  to  yield  a 

maximum  value  of  a  . 

z 

The  various  tests  and  the  specimen  geometries  are  summar¬ 
ized  in  table  8. 

MATERIALS  FABRICATION  AND  QUALITY  CONTROL 

Fabrication  of  the  composite  laminates  was  accomplished 
with  contemporary  vacuum  bag  autoclave  techniques  with  the 
proper  cure  cycle  as  supplied  by  the  material  producer.  In 
order  to  evaluate  process  conditions  for  adequacy,  lamina 
characterization  tests  were  performed  to  determine  fiber  ten¬ 
sion  and  in-plane  shear  properties.  These  properties  are 
summarized  in  table  9.  Review  of  the  lamina  characterization 
tests  indicated  that  the  fabricated  T-300/5208  composite 
possessed  an  average  fiber  tensile  strength  of  1668  MPa,  a 
Young's  modulus  of  141  GPa,  and  an  in-plane  shear  modulus  of 
5  GPa.  These  properties  are  indicative  of  high-quality 
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laminate  fabrication  (ref.  7)  and  suggest  that  the  processing 
conditions  were  adequate. 

Each  composite  panel  was  subjected  to  ultrasonic  inspec¬ 
tion  after  fabrication.  In  this  way,  both  defects  could  be 
isolated  from  test  specimens.  A  typical  "C"  scan  of  a  32-ply 

[±45]  laminate  is  shown  in  figure  31.  The  "C"  scan  reveals 

s  - 

anomolies  near  the  "raw"  panel  edges  associated  with  variations 
in  resin  content  or  thickness  due  to  edge  effects.  These 
regions  of  the  panel  were  removed  prior  to  specimen  fabrica¬ 
tion.  Interior  regions  are  revealed  as  uniform  in  both  thick¬ 
ness  and  resin  content. 

All  panels  were  fabricated  from  a  single  batch  of  the 
T-3 00/52 08  prepreg.  The  panels  were  fabricated  during  the 
initial  phase  of  the  program  and  stored  in  dessicant  for 
future  use.  Specimen  fabrication  was  accomplished  by  first 
cutting  the  composite  panel  into  sublaminates.  The  sublaminate 
was  then  mounted  in  a  fixture  for  bonding  tabs.  After  tabs 
were  bonded  to  the  sublaminates,  test  specimens  were  sawed 
to  the  proper  geometry  with  a  precision  diamond  saw  facility. 
Each  specimen  was  then  engraved  for  identification  as  to 
laminate,  sublaminate,  and  position  within  the  sublaminate. 
Next,  specimens  were  instrumented  with  electrical  resistance 
strain  gages  and  stored  for  subsequent  testing. 

TEST  METHODS 

The  test  methods  developed  in  this  study  include  methods 
for  determination  of  0°  compression  (constrained  and  uncon¬ 
strained)  fatigue  properties,  the  in-plane  shear  fatigue  proper 
ties  in  the  presence  of  a  compressive  transverse  stress,  inter¬ 
laminar  shear  fatigue  properties,  and  interlaminar  normal 
fatigue  properties  of  the  T-300/5208  material  systems. 

Compression  Tests 

The  test  method  chosen  for  evaluation  of  the  0°  compres¬ 
sion  and  in-plane  shear  fatigue  properties  was  a  modified 
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IITRI  wedge  action  grip  test  fixture  shown  in  figure  32.  The 
IITRI  compression  fixture  consists  of  wedge  action  friction 
grips  mounted  in  support  blocks,  with  connecting  rods  and 
linear  bearings  which  yield  a  concentric  load  path.  While 
the  IITRI  fixture  was  originally  designed  for  a  narrow  test 
specimen  of  6.35  mm  in  width  with  an  unsupported  length  of 
6.35  -  12.7  inches,  the  fixture  was  modified  to  accept  speci¬ 
men  widths  up  to  19.05  mm.  In  addition,  specimen  lateral 
deflection  restraint  was  achieved  by  special  roller  supports 
attached  to  the  wedge  action  grips  as  shown  in  figure  33. 

In  the  evaluation  of  0°  compression  strength  of  the 
T-300/5208  composite  system,  the  test  specimen  geometry 
strongly  influences  test  results.  In  addition  to  unsupported 
specimen  length,  the  specimen  width  was  found  to  affect  com¬ 
pressive  strength.  For  example,  the  average  strength  for  a 
specimen  geometry  of  12.7  mm  width  and  28.575  mm  unsupported 
length  was  957  MPa,  while  the  average  strength  for  the  6.35 
mm  width  specimen  was  1134  MPa.  By  reducing  the  unsupported 
length  of  the  6.35  mm  width  specimen  to  15.875  mm, the  average 
compressive  strength  achieved  was  1347  MPa  with  a  corresponding 
modulus  of  133  GPa .  These  results  compared  favorably  with  the 
tensile  properties  (1668  MPa  and  141  GPa)  and,  therefore,  the 
last  specimen  geometry  was  chosen  for  fatigue  testing.  The 
test  specimen's  unsupported  length  is  important  in  the  deter¬ 
mination  of  compressive  strength  due  to  restraint  of  axial 
cracking  rather  than  column  instability  about  the  axis  of  mini¬ 
mum  moment  of  inertia.  Indeed,  the  failure  mode  of  the  0° 
compression  test  specimens  was  axial  splitting  followed  by 
micro-instability  of  the  segments  formed  by  the  axial  cracking. 
It  would  appear,  therefore,  that  the  compressive  properties 
developed  in  this  effort  are  uniquely  a  function  of  the  test 
specimen  geometry.  The  actual  compressive  strength  of  the 
composite  specimen  is  bonded  by  the  Euler  critical  load  for  the 
total  specimen  and  the  strength  achieved  when  axial  cracking 
is  fully  constrained. 
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To  investigate  the  influence  of  axial  cracking  upon  the 

compressive  strength  of  the  composite,  the  [0/90/0] 2g 

laminate  was  found  to  exceed  a  "rule  of  mixtures"  prediction 

of  strength  (935  MPa)  and  thereby  illustrate  the  synergistic 

influence  of  the  90°  constraint  layers.  In  order  to  minimize 

the  influence  of  edge  effects  upon  the  behavior  of  the 

[0/90/0]  test  specimen,  a  specimen  width  of  19.05  mm  was 
z  s 

chosen.  Since  it  was  shown  earlier  that  specimen  width 
influenced  the  compressive  strength  of  the  [ 0 ]  test  speci¬ 
men,  it  should  be  expected  that  the  large  width  would  lead 
to  conservative  estimates  of  the  [0/90/0] 2g  specimen  strength. 
Hence,  the  difference  between  observed  strength  and  "rule  of 
mixtures"  estimate  (7%)  could  have  been  greater  had  the  speci¬ 
men  width  been  smaller  and  the  associated  deleterious  influ¬ 
ence  of  edge  effects  been  absent.  In  order  to  evaluate  the 
influence  of  unsupported  length  upon  compressive  strength  of 
the  [0/90/0] 2g  laminate,  specimen  tests  with  unsupported 
lengths  of  12.7  mm  and  9.525  mm  were  conducted.  The  differ¬ 
ence  in  average  strengths  for  the  two  test  lengths  was  less 
than  one  percent.  It  would  appear,  therefore,  that  the 
lateral  support  system  was  effective  over  such  a  variation 
in  unsupported  length  and  acted  to  insure  compressive  failure. 
It  should  be  noted,  however,  that  the  typical  stress-strain 
response  revealed  some  bad  eccentricity  as  exhibited  by  dif¬ 
ferences  in  apparent  moduli  between  stress-strain  responses 
of  gages  on  opposite  sides  of  the  specimen.  Failure  commonly 
corresponded  to  an  abrupt  increase  in  strain  on  one  side  of 
the  specimen  and  a  corresponding  decrease  on  the  other  side. 
This  behavior  would  tend  to  imply  an  instability  induced 
failure.  However,  the  lateral  constraint  at  the  center  of 
the  specimen  should  preclude  gross  instability.  Hence,  the 
failure  of  the  [0/90/0] 2g  specimen  is  considered  to  be  ini¬ 
tiated  by  a  local  instability  of  the  outer  0°  layers. 

In  order  to  investigate  the  in-plane  shear  fatigue  proper¬ 
ties  of  the  T- 300/5208  material  system  in  the  presence  of 
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a  compressive  transverse  normal  stress,  [±4516]s  laminate 
specimens  were  loaded  in  compression.  Since  the  compression 
response  of  the  [±45]  laminate  is  of  nonlinear  softening 
character,  it  was  necessary  to  employ  a  32-ply  laminate  in 
order  to  preclude  instability  induced  failure.  A  specimen 
width  of  19.05  mm  was  chosen  in  order  to  minimize  the  influ¬ 
ence  of  edge  effects  upon  laminate  response.  By  introducing 
a  teflon  layer  between  specimen  and  tab  adjacent  to  the  test 
section,  the  test  section  could  be  effectively  increased  while 
maintaining  the  lateral  constraint  of  the  tab.  In  order  to 
assess  the  influence  of  the  transverse  normal  stress  upon 
static  in-plane  shear  strength,  it  is  informative  to  compare 
the  tensile  and  compressive  strengths  of  the  [+45]  laminate. 
The  average  static  compressive  strength  observed  was  201  MPa 
with  an  average  modulus  of  18.9  GPa .  The  average  static 
tensile  strength  for  the  identical  specimen  geometry  was  170 
MPa  with  a  modulus  of  18.8  GPa.  The  latter  data  show  good 
agreement  with  the  tensile  response  of  the  8-ply  [  ± 4 5 ] s 
laminate  of  25.4  mm  width  and  152.4  mm  test  section  length. 
Average  tensile  strength  and  modulus  for  this  specimen  geometry 
were  found  to  be  150  MPa  and  19.2  GPa,  respectively.  Hence, 
i-t  would  appear  that  there  is  a  strong  dependence  of  shear 
strength  upon  the  sign  of  the  transverse  normal  stress. 

The  [±45,,]  laminate  specimen  was  subjected  to  compres- 
16  s 

sion  fatigue  employing  the  modified  IITRI  test  fixture  shown 
in  figure  33.  Both  fatigue  life  and  residual  strength  data 
were  developed  with  this  test  method.  Residual  strength  test 
specimens  were  instrumented  with  longitudinal  and  transverse 
strain  gages  for  determination  of  residual  modulus  and 
Poisson's  ratio. 


Tension  Tests 

Tensile  test  methods  were  developed  and  evaluated  for 
determination  of  the  interlaminar  shear  fatigue,  and  inter¬ 
laminar  normal  fatigue  properties.  The  interlaminar  shear 
fatigue  test  specimens  consisted  of  a  [±45/0^1 s  laminate 
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wherein  the  0°  layers  were  discontinuous  at  the  center  of  the 
specimen.  Load  carried  by  the  0°  layers  was  transferred 
through  interlaminar  shear.  Hence,  failure  was  initiated  at 
the  discontinuity  of  the  0°  layers  in  an  interlaminar  shear 
mode.  Specimens  were  subjected  to  a  fatigue  loading  (R  =  0.1, 
30  Hz)  and  the  initiation  of  cracking  detected  visually. 

Static  samples  were  instrumented  with  acoustic  emission  sen¬ 
sors  for  detection  of  interlaminar  failure.  Specimens  were 
not  tested  beyond  the  initiation  of  damage. 

The  interlaminar  normal  stress  fatigue  properties  were 
evaluated  with  the  [±252/90]s  tensile  coupon.  Due  to  the 
large  difference  in  Poisson's  ratios  between  the  90°  and  ±25° 
layers,  failure  is  initiated  at  the  specimen  edge  in  an  inter¬ 
laminar  tensile  mode.  Fatigue  specimens  were  subjected  to 
tensile  fatigue  loadings  of  R  =  0.1  at  a  rate  of  30  Hz.  Both 
acoustic  emission  and  visual  detection  techniques  were  employed 
to  establish  failure  initiation. 

DISCUSSION  OF  RESULTS 

As  noted  above,  the  objective  of  the  uniaxial  compression 

static/fatigue  tests  of  [ ± 4 5 ] s  Graphite/epoxy  laminates  is  to 

determine  the  effect  of  interaction  between  in-plane  shear  and 

transverse  stress  (transverse  tension  plus  in-plane  shear  for 

a  tension  test  of  [ ± 4 5 ]  laminate  vs.  transverse  compression 

s 

and  in-plane  shear  for  the  compression  test).  To  that  end, 

static  tension  and  compression  tests  were  conducted  for  the 

[±451  laminate  and  the  results  are  tabulated  in  table  10. 
s 

It  is  seen  that  the  average  compression  strength  is  higher  than 
the  average  tension  strength,  although  the  moduli  are  almost 
identical.  This  is  an  expected  result  and  is  confirmed  by 
analytical  predictions.  Thus  the  [ ± 4 5 ] ^  laminate  cannot  be 
used  realiably  for  predicting  shear  strength  because  the  mode 
of  transverse  stress  has  a  noticeable  effect  on  it. 

Table  11  and  figure  34  show  the  summary  of  the  T-300/5208 
[±45^ laminate  compression  fatigue  data.  Note  that  there 
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is  only  a  slight  wearout  of  strength  and  stiffness.  This  is 
contrary  to  the  significant  wearout  observed  for  the  tension 
fatigue  tests.  Note  also  that  the  scatter  is  rather  high 
for  the  lifetime  data  at  S  =  0.67.  Figures  35  and  36 
illustrate  the  failure  modes  in  static  tension,  and  static 
compression  after  fatigue,  respectively.  The  failure  in  ten¬ 
sion  occurs  primarily  due  to  ply  separation  which  is  precipi¬ 
tated  by  transverse  tension  failure  in  the  ±45  laminae.  Micro- 
buckling  at  the  lamina  level  is  the  predominant  failure  mode 
for  compression.  An  interesting  failure  mode  is  observed  in 
figure  37  for  a  [±45]  laminate  subjected  to  compression.  The 
surface  ply  has  buckled  in  figure  37  resulting  in  further  loss 
of  stiffness  and  subsequent  failure. 

The  static/fatigue  data  for  interlaminar  normal  stress 

a  as  obtained  from  a  [±25n/90]  laminate  are  shown  in  tables 
z  27  s 

12  and  13  and  in  figure  38.  In  table  12,  the  initiation  of 
failure  at  the  free  edges  was  monitored  by  acoustic  emission 
and  from  the  load-displacement  chart.  In  figure  38,  the 
scatter  is  rather  large  for  S  =  0.6.  However,  this  is  to  be 
expected.  The  slope  of  the  S-N  curve  does  indicate  a  wearout 
rate  that  is  consistent  with  a  matrix-dominated  behavior.  The 
edge  failure  mode  for  this  specimen  is  illustrated  in  figure  39. 

The  static/fatigue  test  data  for  the  [±45/0^]s  partially 
notched  laminates  are  given  in  tables  14  and  15  and  shown  in 
figure  40.  Failure  is  initiated  in  the  interlaminar  shear 
mode  and  this  is  evident  from  the  C-scan  in  figure  41  and  the 
photomicrograph  in  figure  42.  Referring  to  the  derivations 
in  Appendix  B  and  to  table  Bl,  the  analytical  prediction  based 
on  an  interlaminar  shear  strength  of  68.95  MPa  is  substantiated 
by  this  data.  Interestingly  enough,  the  scatter  for  S  =  0.8 
is  rather  significant  while  for  S =  0 . 5 ,  the  data  points  are 
rather  close.  Because  of  the  notch  blunting  effect  during 
fatigue,  the  residual  strength  for  the  specimen  at  S  =  0 . 6 
is  higher  than  the  average  static  strength.  Figures  43  and 
44  illustrate  the  static  and  residual  strength,  and  fatigue 
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failures,  respectively.  The  fatigue  failures,  however,  are 
visible  only  under  magnification. 

The  static  compression  results  for  the  [0]  and  [0/90/0] 
laminates  are  given  in  tables  16  and  17.  The  compression 
strength  for  the  [0]  laminate  is  consistent  with  the  values 
reported  in  the  literature.  The  modulus  is  also  comparable 
to  the  [0]  tension  modulus  as  reported  earlier  in  table  9. 

The  compression  strength  and  modulus  for  the  [0/90/0] 
laminate  is  comparable  to  the  values  obtained  by  multiplying 
the  corresponding  quantities  for  the  [0]  laminate  by  0.66 
(percentage  of  0  degree  laminae  in  the  laminate) . 

Lifetime  and  residual  strength  data  for  the  [0]  laminate 
are  tabulated  in  table  18  and  plotted  in  figure  45.  An  impor¬ 
tant  feature  of  the  [0]  fatigue  tests  is  that  for  most  of  the 
specimens,  failures  occurred  by  the  separation  of  tabs  because 
of  the  high  stress  levels.  Hence  only  two  fatigue  failures 
were  obtained.  When  a  specimen  experienced  separation  of  tabs, 
it  was  removed  from  the  fixture,  the  tabs  were  rebonded,  and 
the  specimen  was  tested  statically  to  failure.  Hence  the 
result  is  a  residual  strength  data  point.  Because  of  these 
difficulties,  an  S-N  curve  for  the  [0]  laminate  could  not  be 
obtained.  Note  also  that  the  residual  strength  actually  shows 
an  increase  in  strength.  This  appears  to  be  paradoxical.  A 
probable  explanation  for  this  inconsistency  is  the  scatter  in 
the  static  strength  data.  The  compression  fatigue  failures 
for  the  [0]  laminate  are  illustrated  in  figure  46.  The  failure 
sequence  is  splitting  in  the  axial  direction,  instability, 
and  brooming. 

The  [0/90/0] 2S  compression  fatigue  test  results  are  shown 

in  table  19  and  figures  47  and  48.  Once  again,  the  data  points 

for  S  =  0.8  and  0.67  sometimes  differ  by  two  or  three  decades. 

While  the  residual  strength  for  the  [0]  shows  an  increase  for 

S  =  0.5,  that  for  the  [0/90/0]  laminates  registers  a  small 

^  s 

decrease.  Note,  however,  that  the  failure  modes  for  [0]  and 
[0/90/0] 2S  laminates  are  different.  In  figure  49,  the  residual 
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strength  test  failure  suggests  the  buckling  of  the  surface  0 
degree  layer  on  an  elastic  foundation  followed  by  overall 
instability  and  subsequent  fracture.  The  fatigue  failures  in 
figure  49  for  S  =  0.8  occurred  generally  inside  the  tabs, 
suggesting  that  some  of  the  data  points  in  figure  47  are  pre¬ 
mature  failures. 


28 


CONCLUDING  REMARKS 


The  basic  failure  model  for  predicting  through-the- 
thickness  crack  growth  and  ultimate  failure  of  notched  com¬ 
posite  laminates  (ref.  3)  has  been  expanded  to  include  the 
effects  of  interlaminar  stresses  in  the  vicinity  of  the  notch. 
The  modification  was  accomplished  by  locally  discretizing  the 
laminate  into  laminae  in  the  region  of  the  projection  of  the 
notch  diameter  along  the  direction  of  the  loading.  In-plane 
and  interlaminar  damage  (inelastic  region  or  crack)  growth 
was  calculated  for  various  laminates  ( [0o/±45]  ,  [0„/±45/0] 

and  [0/±45/90]g)  and  stacking  sequences. 

Because  of  the  possibility  of  a  large  combination  of 
in-plane  and  interlaminar  damage  zones  and  the  numerical  com¬ 
plexity  associated  with  "tracking"  them,  the  capability  of 
the  failure  model  is  restricted  to  the  prediction  of  only  a 
limited  number  of  combinations  of  these  regions.  In  order 
to  simplify  the  analysis,  a  secant  modulus  approximation  to 
the  nonlinear  in-plane  and  interlaminar  shear  stress-strain 
curves  was  considered  in  lieu  of  the  elastic,  perfectly 
plastic  approximation.  This  eliminated  the  need  to  consider 
inelastic  regions  and  damage  was  predicted  as  a  crack  in  a 
layer  or  at  the  interface.  This  enhanced  the  capability  of 
the  analysis  to  predict  axial  crack  growth  in  an  unstable 
fashion  or  until  failure. 

Three-dimensional  finite  element  analysis  was  performed 
in  order  to  develop  an  understanding  of  the  nature  of  the 
stress  distribution  around  a  hole  in  [0/±45]s  and  [±45/0] 
laminates.  Ths  SAP-IV  computer  program  was  utilized  for 
this  purpose.  The  results  obtained  from  the  finite  element 
analysis  were  compared  with  that  of  the  approximate  stress 
analysis  procedure  with  the  objective  of  qualitatively  assess¬ 
ing  the  kinematic  assumptions. 

As  an  additional  task,  a  2-D  elastic-plastic  finite 
element  analysis  was  conducted  for  some  notched  laminates. 
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Because  of  the  consideration  of  inelastic  effects  in  shear  in 
the  basic  failure  model,  it  appeared  desirable  to  conduct  such 
an  analysis  for  the  purpose  of  investigating  the  nature  of  the 
damage  growth  and  the  stress  distribution  in  the  vicinity  of 
the  notch,  and  to  ascertain  the  effect  of  laminate  anisotropy 
on  damage  propagation. 

The  experimental  program  basically  consisted  of  the  fol¬ 
lowing  exploratory  tests: 

(i)  [ ±  4 5 ] s  laminate  uniaxial  compression/compression 

fatigue  tests  to  determine  the  nature  of  the  interac¬ 
tion  between  in-plane  shear  and  transverse  compres¬ 
sion  in  a  lamina  and  to  compare  the  wearout  rates 

with  that  of  a  [±45]  uniaxial  tension/tension 

s 

fatigue  test.  The  results  indicated  that  there  is 
no  noticeable  wearout,  thus  suggesting  that  the 
sign  of  the  transverse  stress  has  an  important  effect. 


(ii)  [0./90.]  laminate  uniaxial  compression/compression 

i  1  s 

fatigue  tests  to  ascertain  the  effect  of  transverse 
restraint  (provided  by  the  90  degree  laminae)  on  the 
failure  mode  and  the  wearout  rate  for  a  unidirec¬ 
tional  lamina  in  compression/compression  fatigue. 

The  data  scatter  for  the  [0/90/0] 2g  laminate 
fatigue  tests  was  rather  large  for  S  =  0.8  and 
specimen  failures  were  experienced  outside  the  gage 
section.  Problems  in  the  [0]  laminate  tests  were 
premature  failures  along  the  bond  between  the  speci¬ 
men  and  the  tabs,  and  an  increase  in  the  residual 
strength  after  fatigue.  The  former  could  be  attrib¬ 
uted  to  high  stress  level  while  the  latter  could  be 
a  result  of  a  scatter  in  the  static  strength  data. 
Because  of  these  inconsistencies,  an  appropriate 
evaluation  of  the  [0/90/0] 2s  laminate  to  obtain 
unidirectional  compression  fatigue  data  could  not 
be  made . 
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(iii)  [±252/90]  laminate  specimen  to  obtain  the  S-N 
curve  for  interlaminar  normal  tensile  stress. 

(iv)  Partially  notched  [ ± 4 5 . / 0  .  ]  laminate  for  deter- 

1  D  s 

mining  the  behavior  of  interlaminar  shear  stress 
in  fatigue. 

The  current  analysis  includes  the  effects  of  damage  growth 
due  to  both  in-plane  and  interlaminar  shear  stress  behavior. 

The  transverse  stresses  in  the  notched  laminate  are  ignored 
in  the  analysis  in  comparison  to  the  axial  stresses.  This 
assumption  was  validated  by  finite  element  analyses  and  experi¬ 
mental  observations.  The  experimental  observations  on  a 
notched  [02/±45]s  laminate  were  predicted  qualitatively  by 
this  analytical  model  using  an  elastic,  secant  modulus  approach. 

The  shortcoming  of  the  failure  prediction  methodology  is 
the  extensive  effort  needed  to  automate  the  numerical  proce¬ 
dure  to  account  for  all  possible  damage  sequences  in  a  general 
laminate.  A  recommended  alternative,  for  future  research, 
would  be  the  utilization  of  these  defined  damage  and  failure 
modes  in  a  finite  element  analysis.  This  could  provide  a 
better  understanding  of  the  stress  state  around  the  notch  and 
have  the  additional  advantage  of  incorporating  all  the  stresses 
for  general  in-plane  loading  into  the  failure  prediction 
methodology . 
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APPENDIX  A 

FAILURE  MODEL  FOR  NOTCHED  COMPOSITE  LAMINATES 
WITH  INTERLAMINAR  EFFECTS 

Referring  to  figures  1  and  2  in  the  main  text,  the  govern¬ 
ing  equilibrium  equations  for  the  various  regions  in  the  notched 
laminate  may  be  written  as  shown  below. 


ELASTIC  REGION 

The  basic  approach  utilized  in  formulating  the  present  model 
is  to  discretize  the  laminate  into  laminae  only  in  the  region  of 
the  projection  of  the  notch  width.  The  remaining  portion  of  the 
homogeneous  laminate  is  divided  into  an  average  stress  concentra¬ 
tion  region  and  a  uniform  stress  region  as  in  reference  2  .  It  is 
also  assumed  that  the  symmetric  laminate  has  a  maximum  of  six 
laminae.  Equilibrium  of  each  ply  in  the  central  core  region  and 
the  adjacent  overstressed  region  yields  the  following  equations: 
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Note  that  the  inplane  shear  stress  transfer  region  is  of 
width  " aQ "  for  all  layers  and  the  interlaminar  shear  stress  trans¬ 
fer  is  assumed  to  take  place  over  a  thickness  "d" .  Axial  force 
equilibrium  in  the  laminate  requires  that: 
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Integrating  equation  (A5)  with  respect  to  x  yields: 
2a2hExu2u2(x)  =  Fx  +  C10  '  2ahlExlu0l(x) 


2ah2Ex2U02(x)  -  ah3Ex3U03(X>  '  2a0hExulVx) 


(A6) 


where  C1Q  is  a  constant  of  integration. 

Substituting  equation  (A6)  into  (A4) ,  the  governing  equations 

can  be  written  in  a  matrix  form  as: 
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and  the  subscript  or  a  superscipt  "e"  represents  the  elastic 
region. 

The  general  solution  of  equation  (A7)  is: 
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where  { S  } e  =  —  [ B ]  e  ^{D},  (AlO) 

m1,m2,m3  and  m4  are  the  eigenvalues  of  the  matrix  [ B] e , 

(u }  (1)  ,  {U}(2),  {U}(3)  and  {U}(4)  are  the  corresponding 

eigenvectors 

and  cio,Cll'C12'  etc*'  are  constants* 

Finiteness  of  the  stresses  for  large  x  is  introduced  by- 
setting  C16,  C 17,  C18  and  C19  equal  to  zero.  Hence 
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Five  boundary  conditions  are  required  to  solve  for  the  five  con¬ 
stants.  If  the  applied  load  is  such  that  the  stresses  in  the 
laminate  are  in  the  elastic  range,  these  conditions  are  the  zero 
stress  condition  on  the  notch  surface  (x  =  0)  and  symmetry 
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conditions  in  the  u^  and  regions  at  the  same  location.  That 
is. 


du0i(0) 
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INELASTIC  REGION (S)  (0<x  <  a) 


Three  inplane  and  two  interlaminar  inelastic  shear  stress 
transfer  regions  are  postulated  in  the  discretized  region.  If, 
for  example,  the  top  layer  (layer  number  1)  has  an  inplane 
inelastic  region  of  length  a  (see  fig.  1  ),  equations  (Al)  and 
(A4 )  are  modified  in  accordance  with  the  elastic-perfectly  plastic 
shear  behavior  assumption  (fig.  6  ) : 
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Equations  (A2) ,  (A3)  and  (A5)  remain  unaltered. 

For  an  inelastic  region  at  an  interface  (see  fig.  1  ) ,  a 
similar  modification  is  introduced.  Assuming  an  interlaminar 
inelastic  region  between  layers  1  and  2,  for  example,  equations 
(Al)  and  (A2)  are  modified  as  shown  below: 
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(A15) 
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Equations  (A3),  (A4)  and  (A5)  remain  the  same. 

The  governing  equations  for  a  region  of  the  laminate  where 
inelasticity  is  present  may  then  be  written  as: 

d2{u} 

- 2^  "  [B]p{u}p  =  (Fx  +  C5^d}  +  (A17) 

dx  F  F 

where  {R}  represents  the  inelastic  contribution  to  the  nonhomo- 
geneous  term.  The  solution  for  this  set  of  equations  is: 
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where  {T}  =  — [B]  1 { R }  and  { S }  =  — [B]  1 { D }  (A19) 

P  P  P 

and  the  subscript  or  the  superscript  "p"  represents  the  inelastic 
region . 

The  eigenvalues  m.(i  =  5,8)  and  the  corresponding  eigenvectors 

.V  1 

{U}  of  [B]  for  the  inelastic  region  are  different  from  those 
P 

for  the  elastic  region  due  to  the  changes  in  the  [B]  matrix. 
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Coordinate  x  for  each  region  is  measured  from  the  origin  of  that 
region. 

Consider  now  an  inelastic  region  in  the  laminate,  extending 
over  a  length  of  a  from  the  notch.  Equations  (All)  and  (A18) 
are  the  displacement  solutions  for  the  elastic  and  inelastic 
regions.  The  boundary  conditions  for  the  problem  are: 
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In  addition  to  the  above  14  boundary  conditions,  the  inelas¬ 
ticity  condition  for  the  appropriate  layer/interface  at  the  junc¬ 
tion  of  the  two  regions  must  be  imposed  for  a  unique  F  -  a 
relationship . 

REGION  WITH  AN  INPLANE  OR  AN  INTERLAMINAR  CRACK  (0  <  x  <  C) 

The  governing  equations  for  a  region  with  a  crack  are 
obtained  by  imposing  the  stress-free  conditions  on  the  cracked 
inplane  or  interlaminar  surface.  For  an  inplane  crack  in  the 
top  layer,  for  example,  equations  (Al)  and  (A4)  are  modified  to: 
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(A21) 
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Equations  (A2),  (A3)  and  (A5)  remain  unchanged. 

If,  on  the  other  hand,  there  is  a  debond  between  layers  1 
and  2,  equations  (Al)  and  (A2)  take  the  form: 
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Equations  (A3) ,  (A4)  and  (A5)  remain  unchanged. 

For  either  type  of  crack,  the  governing  equations  take  the 
form  of  equation  (A7) ,  and  the  solutions,  the  form  of  equation 
(A9) .  Introducing  new  constants, 
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where  (S)  =  — [B]  ^  D  (A27) 

c  c 

and  the  superscript  or  subscript  "c"  represents  the  cracked 
region . 

The  boundary  conditions  for  a  problem  with  a  cracked  region, 
an  inelastic  region,  and  an  elastic  region  could  be  written  in  a 
manner  similar  to  equations  (A20) .  For  a  unique  F  —  a  -  ? 
relationship,  the  criterion  for  crack  initiation  will  also  have 

to  be  imposed  at  the  crack  tip. 

DEGENERATE  CASES 

The  various  regions  in  the  laminate  as  described  above  are 
determined  by  the  sequence  in  which  defects  are  likely  to  be 
initiated  with  an  increase  in  the  applied  load.  However,  more 
than  one  defect  could  initiate  and  grow  simultaneously.  For 
some  widely  used  laminate  constructions,  certain  sequences  of 
defect  formation  could  cause  the  nature  of  the  above  solutions 
to  change.  Such  a  situation  will  arise  (i)  when  two  or  more 
inelastic  zones  ( inplane/interlaminar )  appear  simultaneously,  or 
(ii)  when  one  or  more  laminae  act  as  non-load  carrying  members 
undergoing  rigid  body  displacements.  The  first  type  of  degeneracy 
occurs  for  elastic-plastic  behavior  in  inplane/interlaminar  shear, 
and  the  degenerate  layer/layers  have  a  linear  variation  of  axial 
strain.  The  second  type  of  degeneracy  is  more  likely  to  occur 
in  an  elastic  analysis  with  a  shear  behavior  governed  by  a  secant 
modulus  approach.  The  degenerate  layer/layers  in  this  case  are 
essentially  "carried"  by  the  load-carrying  adjacent  regions. 
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Simultaneous  Inelastic  Inplane  and  Interlaminar  Shear  Stress 
-  —  Transfer  for a  Lamina 

Consider,  for  example,  the  first  type  of  degeneracy  when 
inelastic  zones  have  initiated  simultaneously  in  the  surface 
layer  1  (inplane  shear  mode)  and  at  the  interface  between  layers 
1  and  2  (interlaminar  shear  mode).  Equation  (Al)  takes  the  form: 
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Integrating  this  equation  with  respect  to  x  twice,  one  obtains 
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If  the  boundary  of  this  region  for  x  =  0  is  the  notch  surface, 
then  C^q  vanishes.  Equation  (A 2 9)  then  becomes: 
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(A3  0 ) 


Equation  (A2 )  is  then  modified  to: 
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and  equation  (A4)  assumes  the  form: 
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Upon  utilizing  the  solution  in  equation  (A30) ,  equation  (A6) 
becomes 
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Substituting  equation  (A33)  into  equation  (A32) ,  and  considering 
equations  (A31)  and  (A3),  the  governing  equations  for  the 
degenerate  region  are  written  as: 
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and  the  subscript  or  superscript  "d"  refers  to  the  degenerate 
solution(s).  The  solution  to  equation  (A34)  is: 
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The  boundary  conditions  for  a  laminate  with  a  degenerate 
region  are  written  in  the  same  form  as  before. 

Simultaneous  Inplane  and  Interlaminar  Cracking 

for  a  Lamina 

If,  instead  of  the  elastic-plastic  approach,  an  elastic 
secant  modulus  approach  is  chosen  to  predict  the  "peeling  off" 
of  the  surface  lamina (e)  in  a  laminate,  the  above  solutions  could 
be  modified  accordingly.  Considering  the  same  example  presented 
earlier,  the  surface  ply  degeneracy  of  the  second  type  changes 
equations  (A30)  through  (A33)  by  eliminating  all  the  terms  with 
yield  stresses.  The  governing  equations  then  reduce  to  the 
form: 
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The  solution  to  equation  (A40)  is: 
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It  is  seen  that  a  single  ply  degeneracy  reduces  the  number 

of  coupled  equations  by  one,  as  reflected  by  the  matrix  [ B ] d . 

If  a  group  of  layers  has  a  degenerate  solution  in  a  region  of 

the  laminate,  the  [B]  ,  matrix  reduces  in  size  further.  The  equa- 

d 

tions  and  the  solutions  for  these  cases  could  then  be  written 
as  above. 


Simultaneous  Inplane  Plasticity  in  All  Layers 


An  extreme  case  of  inplane  inelasticity  in  all  the  layers 
over  a  region  (0  <  x  <  a)  of  the  laminate  is  not  uncommon.  In 
that  region,  the  core  region  (see  fig.  1  )  experiences  a  plastic 

flow.  Mathematically,  this  uncouples  the  core  displacement  equa¬ 
tions  from  those  for  the  u^  and  u^  displacements.  Equation  (A4) 
for  this  region  is: 
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where  b3  =  2x^1^  +  2ixy2h2  +  ixy3h3 


(A4  3 ) 


Equation  (A6)  becomes: 
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The  last  two  equations  can  be  solved  simultaneously  for  u^ (x) 
and  u2 (x)  in  the  degenerate  region  of  the  laminate. 

The  core  displacements  in  this  region  are  governed  by  the 
plastic  flow  equations.  Each  layer  has  an  inplane  shear  force 
that  is  proportional  to  a,  and  an  interlaminar  shear  strain  that 
maintains  the  value  at  x  =  a.  Therefore,  the  total  axial  force 
in  each  layer,  or  the  axial  strain  in  that  layer,  varies 
linearly  with  respect  to  x.  That  is, 
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x2  22  xy2  2 


Txzl2(a)a  +  TXZ23<a)a 


(A47) 


and  E  ah_E  =  2x  h_ 
x3  33  xy3  3 


2  x  _  ( o, )  a 

xz23 


(A48) 


The  core  displacement  solutions  for  the  degenerate  region 
are  dependent  on  the  displacement  solutions  of  the  adjacent  region 

(along  the  x-axis)  through  the  TXZ]_2^a)  and  Txz23^  terms- 
In  other  words, 
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x  „ (a) /degenerate  region  =  x  , 9 ( 0 ) /adjacent  region 
xzl/  xzxz 

x  ( a) /degenerate  region  =  x  (0)/adjacent  region 

XZ23  XZz j 


(A49) 


When  an  elastic  analysis  with  a  secant  modulus  shear  behavior 
is  carried  out,  a  region  (0  <  x  <  ?)  of  through-the-thickness 
crack  might  be  predicted.  For  this  region,  equations  (A42)  and 
(A44)  become: 


(u.  -  u2) 

Exula? 


(A50) 


2a„hE  „u„ 
2  xu2  2 


Fx  +  C 


-  2a  hE  . u,  , 
0  xul  1 


(A51) 


respectively.  The  above  equations  yield  solutions  for  u^ (x)  and 
U2 (x)  for  the  region  0  <  x  <  C. 

The  core,  in  this  case,  moves  as  a  rigid  body,  and  the  dis¬ 
placements  are  those  of  the  adjacent  region  along  the  x-axis. 

UQi/degenerate  region  =  uQi ( 0 ) /ad jacent  region,  i  =  1,3  (A52) 
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APPENDIX  B 
TEST  SPECIMEN  DESIGN 


The  objective  of  the  experimental  program  is  primarily  to 
perform  exploratory  tests  for  determining  basic  data  required  as 
input  into  the  failure  model.  After  considering  the  deficiencies 
in  the  data  during  the  previous  analytical/experimental  correla¬ 
tion  study  (ref.  2  )  and  the  needs  for  the  modified  basic  failure 

model  with  interlaminar  effects,  the  following  tests  were  identi¬ 
fied  : 

(i)  [±45]  laminate  uniaxial  compression/compression  fatigue 

tests  to  determine  the  nature  of  the  interaction  between  inplane 
shear  and  transverse  compression  in  a  lamina  and  to  compare  the 
wearout  rate  with  that  of  a  [±45]  uniaxial  tension/ tension 
fatigue  test  ; 

(ii)  [ 0 . / 9 0  .  ]  laminate  unaxial  compression/compression 
i  1  s 

fatigue  tests  to  ascertain  the  effect  of  transverse  restraint 
(provided  by  the  90  degree  laminae)  on  the  failure  mode  and  the 
wearout  rate  for  a  unidirectional  lamina  in  compression/compres¬ 
sion  fatigue  ; 

(iii)  [±25„/90]  laminate  specimen  to  obtain  the  S-N  curve 
2  s 

for  interlaminar  normal  tensile  stress  and 

(iv)  Partially  notched  [ ± 4 5 . / 0  .  ]  laminate  for  determining 

i  1  s 

the  behavior  of  interlaminar  shear  stress  in  fatigue. 

With  regard  to  the  design  of  all  specimens  subjected  to  com¬ 
pression,  dimensions  were  chosen  such  that  failure  would  not 

occur  prematurely  due  to  instability.  For  the  [0./90-]  speci- 

—  j  s 

mens,  the  stacking  sequence  was  [0/90/0] ^  .  This  would  avoid 

any  interlaminar  edge  failures  due  to  compression  loading.  The 

design  of  the  [±25/90]  laminate  was  taken  from  reference 

2  s 

This  particular  stacking  sequence  results  in  a  high  value  of  the 
interlaminar  normal  tensile  stress  when  the  specimen  is  subjected 
to  tension  loading.  This  ensures  edge  failure  due  to  interlam¬ 
inar  normal  stress.  Finally,  a  partially  notched  [ ± 4 5 . / 0 . ] 

i  j  s 
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laminate,  with  the  notch  in  the  0  degree  layers,  was  designed 
to  investigate  the  interlaminar  shear  behavior. 

Figure  B1  illustrates  a  partially  notched  [±45^/0j]g  lamin¬ 
ate  with  the  notch  in  the  0  degree  layers  across  the  width  of 
the  specimen.  Consideration  of  equilibrium  in  the  x-direction 
results  in  the  following  differential  equations: 


2 

E1h1  d  ux  u2~ui 

~2  +  ~ d 

dx 


xz  =  0 


(Bl) 


E„h0  d  u2 

2  2  ,  2 

dx 


U2_U1  G 


xz  =  0 


(B2) 


where  'd'  is  the  width  of  the  interlaminar  shear  transfer  zone. 
Subscripts  1  and  2  correspond  to  the  0  degree  and  ±45  degree 
plies,  respectively.  Also,  force  equilibrium  at  any  cross  sec¬ 
tion  of  the  laminate  results  in  the  following  expression  for 
the  displacement  u2 : 


rFx  „ 

u„  =  [tt—  -  C,  - 
2  2w  1 


ElhlUl- 


E2h2 


(B3) 


ElhlUl (0) 

where  the  constant  *  -  - 2 -  and  F  as  the  applied 

laminate  force.  Upon  substitution  of  equation  (B3)  in  equa¬ 
tion  (Bl)  and  the  subsequent  solution  of  the  resulting  dif¬ 
ferential  equation,  one  obtains 


—  ry  v 

u1  =  C2e  +  Ax  +  B 


(B4) 
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where  is  a  constant  of  integration. 


2 

a 


2G 


xz 


Elhld 


Eh 

[1  +  2ihH 

Zh2n2 


(B5) 


FG 

A  = - — - ~ -  and  B 

WE1E2h  h  a'  d 


-  2C  G 
_ 1  xz 

E, E„h. h_da 


2 


(B6 ) 


The  boundary  condition  ensuring  the  finiteness  of  the  force  for 
large  x  has  already  been  considered  by  eliminating  exponential 

(XX 

terms  of  the  type  e  .  Consideration  of  the  boundary  condition 
du 

gj  =  0  at  x  =  0  yields 

C2  =  A/a  (B7 ) 

and  consequently, 

-FG 

C,  =  - — - 2 -  (B8) 

2Wa (Enh„da  -G  ) 

2  2  x  z 


An  expression  for  the  displacement  u2  can  be  obtained  from 
equations  (B4) ,  (B6)  and  (B3) .  Having  known  the  displacements 

u^  and  u2 ,  the  interlaminar  shear  stress  can  be  calculated  as 


L 

xz 


(X) 


G 

xz 


u  (x) -u 
[— - -r- 


(B9 ) 


The  objective  of  the  above  analysis  is  to  determine  the 
percentage  of  the  0  degree  layers  in  the  [±45^/0^]^  laminate 
such  that  failure  occurs  due  to  the  interlaminar  stress  t 


xz 


in  the  vicinity  of  the  notch  and  not  due  to  overstress  in 
the  ±45°  layers .  Also  failure  should  not  occur  in  the 
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unnotched  portion  of  the  laminate. 

Table  (Bl)  shows  the  results  of  the  analysis.  Assuming 

that  the  ultimate  value  of  the  interlaminar  stress  x  is 

2  XZ 
68.94  MN/M  (10  ksi) ,  the  40%/60%  [±45/0]s  laminate  appears 

to  be  the  logical  choice.  A  problem  associated  with  the  stack¬ 
ing  sequence  for  this  laminate  is  that  for  tension  loading,  the 
interlaminar  normal  stress  at  the  free  edge  would  be  tensile. 
The  effect  of  a  tensile  a  can  be  alleviated  by  choosing  a 
wider  specimen  (say,  3.81  cm  or  1-1/2  in)  . 
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Table  1. 


Unidirectional  T-300/5208  Lamina 
Properties  and  Some  Laminate  Data 
Used  in  the  Current  Study 


Elastic-Plastic 

Approach 

Secant  Modulus 
Approach 

Property 

Trial  Data 

Actual  Data 

Actual  Data 

G  (0°)  in  GPa 
xy 

6.550 

6.550 

3.283 

G  (±45°)  in  GPa 
xy 

35.935 

35.935 

35.935 

G  (90°)  in  GPa 
xy 

6.550 

6.550 

3.283 

G  in  GPa 

xy 

3.275 

5.516 

5.516 

Ex(0°)  in  GPa 

137.896 

137.896 

137.896 

Ex (±45°)  in  GPa 

22.394 

22.394 

22.394 

Ex(90°)  in  GPa 

10.342 

10.342 

10.342 

a  (mm) 

6.350 

6.350 

6.350 

(mm) 

2.540 

3.332 

3.332 

W  (mm) 

40.640 

38.100 

38.100 

d=d^=d2  (mm) 

0.070 

0.014 

0.014 

x  __  in  MPa 
xz 

68.948 

68.948 

68.948 

t  '  (0°)  in  MPa 

68.948 

68.948 

68.948 

if  (±45°)  in  MPa 

xy 

448.162 

723.954 

723.954 

_  <0 

Txy  (90°)  in  MPa 

68.948 

68.948 

873.709 

am  ( 0/±45/90)  in  MPa 

4481.620 

482.636 

482.636 

oT(0/0/±45)s  in  MPa 

4481.620 

689.480 

689.480 

YU  (0°) 

xy 

0.021 

0.021 

- 

YUV  (±45°) 
xy 

0.1 

0.0206 

- 

Yu  (90°) 
xy 

1.29 

1.6914 

- 

U 

Yxz 

0.025 

0.025 

- 
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Table  2.  Stresses  in  a  [02/±45/0]s  T-300/5208  Laminate 

For  Different  Interlaminar  Boundary  Layer  Thicknesses 


d  (mm) 

0.070  ' 

0.035 

0.014 

axyl, 3>0 

Txyl  (MPa) 

68.95 

68 . 95 

68.95 

Txy2  (MPa) 

366.78 

372.21 

375.57 

Txy3  (MPa) 

67.51 

68.18 

68.60 

Txzl2  (MPa) 

19.15 

20.22 

21.17 

Txz2  3  (MPa) 

-  5.97 

-5.98 

-5.93 

SCF  (Stress 

Concentration 

Factor) 

2.113 

2.116 

2.117 

a  ,  >  0 ;  a  ~  >  0 

xyl , 3  xy2— 

Txyl  (MPa) 

68.95 

68 . 95 

68.95 

Txy2  (MPa) 

448.16 

448.16 

448.16 

Txy3  (MPa) 

68.95 

68 . 95 

68.95 

txz12  (MPa) 

23.48 

24.44 

25.34 

T-xz2  3  (MPa) 

-7.32 

-7.23 

-7.11 

axyl ,  3  (inrn) 

0.52 

0.48 

0.45 

SCF 

2.1098 

2.113 

2.114 

a  0  >  0 ;  a  ~  >  0 ; 

xyl, 3  xy2 

xyl , 3- 

axyl ,  3  (ltlIn> 

1.93 

1.92 

1.92 

axy2  (n”) 

1.30 

1.36 

1.40 

Txyl  IMPa) 

68. 95 

68.95 

68.95 

T  _  (MPa) 
xy2 

448 . 16 

448.16 

448.16 

T  _  (MPa) 

xy3 

68.95 

68.95 

68.95 

Txzl2(HPa) 

36.13 

38.07 

40.11 

txz23  <HPa) 

-11.25 

-11.27 

-11.40 
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Table  3. 


Stresses  in  a  [±45/0_/0]  T-300/5208  Laminate  For 

2.  s 

Different  Interlaminar  Boundary  Layer  Thicknesses 


d  (mm) 


0.07 

0.035 

0.014 

axy2 , 3>  0 

T^llMPa) 

361.07 

369.17 

374.45 

x  (MPa) 

xy2 

68.44 

68.67 

68.84 

T  , (MPa) 
xy3 

68.95 

68.95 

68.95 

Txzl2(MPa) 

-23.99 

-25.41 

-26.63 

t  _  (MPa) 
xz23 

-4.45 

-4.81 

-5.21 

SCF  (Stress 

Concentration 

Factor) 

2.109 

2.114 

2.117 

axy2,3>0;  axy2- 

>0 

Txyi(MPa) 

448.16 

448.16 

448.16 

Txy2(®a) 

68.95 

68.95 

68.95 

t  -(MPa) 
xy3 

68.95 

68.95 

68.95 

Txz12 (MPal 

-29.88 

-30.96 

-32.00 

Txz23lMPa> 

-5.54 

-5.8  6 

-6.18 

a  (mm) 

xy2 , 3 

0.54 

0.49 

0.45 

SCF 

2.106 

2.111 

2.114 

axy2,3>0;  axy2 

>0'^xyl,3— 0 

axy2 , 3  ^nun) 

1.91 

1.91 

1.91 

a  , (mm) 
xyl 

1.24 

1.32 

1.38 

Txyl(MPal 

448.16 

448.16 

448.16 

Txy2 (MPa) 

68.95 

68.95 

68.95 

Txy3(MPa| 

68.95 

68.95 

68.95 

Txzl2<MPa) 

-45.32 

-47.91 

-50.29 

Txz23IMPa) 

-8.41 

-9.05 

-9.60 
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NOTES  FOR  TABLE  2  AND  TABLE  3 


TABLE  2: 

a  ,  width  of  the  average  stress  concentration  region 
=  1.27  mm. 

In-plane  shear  yield  stresses  =  68.95  MPa  for  a  0®  layer 

448.16  MPa  for  a  ±45°  layer 
Interlaminar  shear  yield  stresses  =  68.95  MPa 

Subscripts  1,  2  and  3  refer  to  the  0^,  ±45  and  0  laminae, 
respectively 

.  Stress  Concentration  Factor  (SCF)  is  the  normal  stress 
in  the  average  stress  concentration  region  divided  by 
the  gross  applied  stress 


TABLE  3: 

Subscripts  1,  2  and  3  refer  to  the  ±45,  0 ^  and  (T  laminae, 

respectively . 
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Table  4.  Damage  Growth  in  a  [0/±45/90]s  T-300/5208 

Laminate  for  Different  Interlaminar  Boundary 
Layer  Thicknesses  (Secant  Modulus  Approach) * 


d  =  0. 

00140  cm. 

d  =  0. 

00279  cm. 

d  =  0. 

00699  cm. 

Force 

0 

Force 

C 

Force 

C 

(KN) 

(cm.  ) 

(KN) 

(cm.  ) 

(KN) 

(cm.  ) 

25.94 

0 . 0** 

25.96 

0.0** 

26.03 

0.0** 

27.25 

.0254 

27.23 

.0254 

27.19 

.  0254 

27.60 

.0635 

27.90 

.1016 

27 . 87 

.1016 

28.1-2 

.  1270 

28 . 69 

.2032 

28.67 

.  2032 

29.04 

.2540 

29.41 

.3048 

29.39 

.3048 

29.88 

.  3810 

30.77 

.  5080 

30.78 

.  5080 

*The  applied  laminate  stress,  a^,  at  which  transverse  crack 
propates  from  the  notch  tip  is  assumed  to  be  very  large. 

**Inplane  crack  initiates  in  the  0°  and  ±45°  layers. 
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Table  5.  Input  Data  for  Different  Laminates  for  a  Two-Dimensiona 
Elastic-Plastic  Analysis  (PLANE  Finite  Element  Program) 


Table  6.  Load  at  Initiation  of  Plasticity 
for  Various  Laminates 


P 

cr 


Laminate 


Table  7  . 


Variation  of  Stress  Concentration 
Factor  for  Various  Laminates  as  a 
Function  of  Applied  Load 


Laminate 

Applied  Load 

Stress  in  Element  12 

Applied  Load 

[0] 

20.82 

1.87 

29.15 

1.92 

33 . 32 

2.0 

[0/90] s 

17.54 

1.77 

21.04 

1.77 

28.06 

1 .65 

t°2/±45Is 

30.85 

2.98 

43.19 

2.43 

49 . 37 

2.14 

58 


Table  8.  Various  Tests  and  Specimen  Geometries 


T5 

CD  £  m 

an  o-p 

•H  fC3  *H  CCi 
fO  *H  CD  Q 
M  X  CO 
-P  (0  0)  «J 
BH  H  C 
C  C  Oj-H 
O  D  a  £ 
o  o  m 
£  uc 
D 


£ 

i— I  o 
m  *h  * 

•HO)* 

x  co  ~ 

m  (D  Pli 

£  acn 

D  g  — 

o 

u 


CD  H  O  -P 

£  m  -h  m 

•H  -H  CO  Q 

m  x  co 
u  m  a)  m 

-P  *h  p  £ 

m  c  a-H 

c  d  a  a 


p  m 

m  -p 

£  m 

•Hi — IQ 

a  m 
maw 

rH  P  CO 
POO) 
(D  £  P 
-P  -P 

£  C/3 

H 


CO 

CO 

CD 

£ 

O 

•H 

r£  CD 
Eh  P 
cj> 
X  *H 
-P 
p  m 
-P  Cn 
TS  I 

■H  (j-i 

£ 


**Unsupported  Length  of  Specimen 


Table  9.  NARMCO  T-300/5208  Graphite  Epoxy- 
Unidirectional  Laminate  Tension 


Properties 


Specimen 

E  (GPa) 

X 

V 

xy 

ot (MPa) 

t 

£ 

X 

L2 

137.9 

.293 

1606 . 7 

11,000 

R3 

137.9 

.  313 

1661.1 

11,700 

Ll 

143.4 

.  2.37 

1669.2 

10,930 

R4 

146.9 

.  310 

1735 . 3 

11,220 

Average 

141.6 

.288 

1668.1 

11,210 
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Table  10.  [±45..,]  T-300/5208  Laminate  Static  Test 

16  s 

Data  for  Uniaxial  Tension  and  Compression 


Specimen 

Type  of  Test 

a  (MPa) 

X 

Ex  (GPa) 

V 

xy 

1BR1 

Uniaxial 

Compression 

202.8 

- 

- 

1BR5 

II 

206.3 

19.8 

.7645 

1BR2 

II 

192.5 

18.0 

.7565 

1CC1 

ll 

203.5 

- 

- 

1CC2 

II 

202.0 

- 

1CC3 

II 

202.9 

_ 

— 

Average  201.7 

18.9 

.7605 

1BR6 

Uniaxial 

Tension 

171.7 

19.4 

0  .6795 

6AR1 

ll 

155.1 

19.0 

.8000 

6AR2 

ll 

150.3 

19.5 

.7400 

6AR3 


I  143.8 
Average  155.2 


18.4 

19.0 


.'.7500 

.7420 


Table  11.  T-300/5208  [±45.,]  Laminate  Fatigue 

lbs 

Compression  Data  (R  =  0.1,  f  =  30  Hz) 


' 

Residua! 

Compress 

ion  Test 

Specimen 

S 

N  (cycl 

es) 

oc  (MP< 

X 

a) 

E  .(GPa) 

X 

V 

xy 

1BR3 

.  80 

12, 

490 

1BL1 

.80 

4, 

450 

1BL3 

.80 

7, 

200 

Avg . 

1AC3 

.  67 

2,146, 

700 

1BR4 

.67 

15, 

450 

1CL3 

.  67 

55, 

650 

1AC4 

.  67 

1,308, 

130 

1CL1 

.  67 

500  , 

000 

199  . 

9 

17.58 

.781 

1BL6 

.67 

500, 

000 

195. 

4 

16.55 

.808 

1AL5 

.  67 

500, 

000 

201 . 

7 

17.44 

.  866 

Avg . 

199. 

0 

17.17 

.818 

1CL4 

.67 

150, 

000 

200. 

4 

-- 

— 

1AC5 

.  67 

150, 

000 

197. 

7 

17.65 

.788 

1BL4 

.  67 

150, 

000 

192  . 

_8 

14 . 55 

.827 

Avg . 

197. 

0 

16.13 

.808 

1CL6  - 

.  50 

10,000, 

000 

201. 

5 

18.68 

.819 

1AL6 

.  50 

10,000 , 

000 

200. 

5 

17 . 65 

.871 

1AL4 

.  50 

10,000, 

000 

198  . 

1_ 

17.51 

.777 

Avg . 

200. 

2 

17. 93 

.  822 

1AL2 

.  50 

1,000, 

000 

198. 

0 

17.10 

.827 

1AC6 

.  50 

1,000, 

000 

198  . 

3 

20.20 

.775 

1BL5 

.  50 

1,000, 

000 

194  . 

_2 

18.00 

TO  A 

.700 

n  c  1 

L 

Avg . 

196  . 

8 

18.41 

.  /  D  / 
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Table  12  .  Static  Test  Data  for 

[±25o/90]  T-300/5208 

2  S 

Laminates 


Specimen 


2AR1 


2BL3 


2AL3 


Type  of 
Test 

X  ' 

(MPa) 

t** 

0  , 

X 

(MPa) 

Uniaxial 

206 

179 

Tension 

n 

224 

226 

ii 

228 

210 

Avg.  219 

205 

*Determined  from  Acoustic  Emission 
**Obtained  from  Load-Displacement  Chart 


Table  13  . 


Tension/Tension  Lifetime  Fatigue 
Test  Data  for  [ ± 2 5 9 0 ] ^  T-300/5208 
Laminates  (f  =  30Hz,  R  =  0.1) 


Specimen 

s 

N  (cycles) 

2BL2 

0 . 8 

4510 

2AR2 

0 . 8 

5550 

2AL1 

0 . 8 

3120 

2AL2 

0.7 

7600 

2BR3 

0 . 7 

8400 

2BL1 

0 . 7 

8000 

2AR3 

0 . 6 

48,660 

2BR2 

0  .  6 

205,550 

2BR1 

0 .6 

794,050 

2CR1 

0.6 

2,061,580 

2DL1 

0.6 

155,600 
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Table  14  .  Static  Test  Data  for 

[±45/0  ]  Partially  Notched 
j  s 

T-300/5208  Laminate 


a 


Specimen 


3AR2 


3AL2 


Type  of  Test 


Uniaxial 

Tension 


fl 


1  (MPa) 


104.1 


77.2 


3AL3 


92.4 

Avg.  91.01 


Table  15.  Tension/Tension  Lifetime  Fatigue 

Data  for  [±45/03]  Partially  Notched 
T-300/5208  Laminate  (R  =  0.1,  f  =  30Hz) 


Specimen 

S 

N  (cycles) 

ot  (MPa) 

X 

3AR1 

0 . 8 

530 , 450 

3AL1 

0 . 8 

82 , 120 

3C2 

0 . 8 

8 ,910 

3BL3 

0.7 

504 ,210 

3BL2 

0.7 

2,957,100 

3BR3 

0.7 

617 , 740 

3BR2 

0.6 

4,397,520 

3BL1 

0.6 

5xl06  (runout) 

113* 

3C1 

. . 

0.6 

6 

5x10  (runout) 

- 

*Residual 
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Table  16.  Uniaxial  Compression  Static 
Test  Data  for  [0]  T-300/5208 
Laminate 


a 


Specimen 


7D3 

7D4 

7D5 

7D6 

7E1 

7D7 

7E2 


t 

x 


(MPa) 


E 


1262.4 
1108 . 7 
1234.9 
1421.0 

1538.2 

1367.2 

1611 . 3 


(GPa) 

x 


131.7 

134.4 

135.1 


7E3 


Avg . 


1234.2 

1347.2 


131.7 

133.1 


Table  17 .  Uniaxial  Compression  Static 
Test  Results  for  [0/90/0] g 
T-30 0/520 8  Laminate 


Specimen 

oC  (MPa) 

X 

E  (GPa ) 

X 

V 

xy 

5BC1 

1061 . 5 

87.6 

0.066 

5AC5 

1038.4 

97.9 

0.091 

5AL5 

867.0 

91.0 

0.056 

5AC  4 

1000.8 

93.1 

0.105 

5BL1 

913.2 

91 . 7 

0.067 

5AC6 

1172.1 

_ 

- 

1008.7 

'  j 

92.4 

0.077 
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Table  18.  Uniaxial  Compression/Compression 

Lifetime  and  Residual  Strength/Stiffness 
Data  for  [0]  T-300/5208  Laminate 
(f  =  30Hz,  R  =  10) 


Residual 

Strength 

0C  (MPa) 

X 

Residual 
Modulus 
E-x  (GPa) 

Specimen 

S 

N,  Cycles 

7D14 

o 

00 

1,420 

1418* 

7D15 

•I 

1,470 

FF 

7.D16 

II 

850 

1146* 

7D25 

It 

4,240 

1217* 

7D26 

tl 

4,110 

1656* 

7D27 

tl 

3,350 

1628* 

7D11 

0.67 

io6 

1274* 

7D12 

tl 

94,510 

1243* 

7D13 

it 

4,677,110 

1192* 

7D22 

It 

10,300 

1772* 

7D23 

If 

10,240 

1487* 

7D24 

11 

7,780 

1573* 

7D21 

0.5 

108,200 

FF 

7D8 

ll 

5  x  106 

1288 

- 

7D9 

It 

5  x  106 

1513 

140 

7D10 

It 

5  x  10^ 

1578 

137 

Avg.  1460 

138.5 

7D18 

tl 

106 

1682 

125 

7D19 

II 

106 

1650 

130 

7D20 

II 

106 

1378 

129 

Avg-  1570 

12B 

FF  Fatigue  Failure 

♦Residual  Strength  Tests  Conducted  After  N  Cycles 


69 


Table  19  .  Uniaxial  Compression/Compression  Lifetime 
and  Residual  Strength/Stiffness  Data  for 
[ 0/90/0 ] 2s  T-300/5208  Laminate  (f=30Hz,  R=10) 


Residual 

Strength 

Residual 

Modulus 

E  ( GPa ) 

X 

Poisson  *  s 
Ratio 

V 

xy 

Specimen 

S 

N , Cycles 

(MPa) 

5BC5 

0 . 8 

150 

5BC  4 

n 

640 

5AL1 

ti 

15,800 

5AL4 

ii 

10 

5BC6 

ii 

2,770 

5AL2 

ii 

1,110 

5BL4 

0.67 

5,920 

5BC2 

n 

25,290 

6* 

5AL6 

ii 

5x10 

816 

97 

0.172 

5CL1 

M 

2,009,050 

1061 

79 

0.05 

5BL2 

ii 

5x10  6 

966 

93 

0.03 

5AL3 

ii 

2,470 

- 

- 

Avg.  1014 

86 

0.04 

6  * 

5BL6 

0.5 

5x10 

829 

83 

0 .183 

5AC3 

ii 

6* 

5x10 

872 

87 

0.06 

6* 

5BL3 

n 

5x10 

988 

8J7 

0.092 

Avg.  897 

86 

0.112 

5AL1 

ii 

io6 

852 

97 

0.072 

5BL5 

ii 

106 

865 

93 

0.044 

5AC2 

H 

106 

842 

9_0 

0.109 

*Runout 

853 

93 

0.075 
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Table  Bl .  Design  of  a  Partially  Notched  [ ± 45 . /0 • ] s 
Laminate  to  Assure  Failure  due  to  1  ^ 

Interlaminar  Stress  t 

X  z 


%  of  0°  Layers 

40% 

60% 

80% 

Initial  Failure  Stress 
for  Laminate  (MPa) 

(±45°  Layers  in  Shear) 

393.0 

554.7 

730.8 

Ultimate  Failure  Stress 
for  Laminate  from 

Netting  Analysis  (MPa) 

(0°  Layer  in  Tension) 

579.2 

868.7 

1158.3 

Failure  Stress  of  ±45° 
Laminate  in  the  Notched 

Area  (MPa) 

124.1 

124.1 

124.1 

Corresponding  Force 
^±45° 

(Specimen  Width=25 . 4mm) 

2642.2 

1761.5 

880.7 

Value  of  the  Force,  Fg, 

when  Failure  Occurs  by 

Interlaminar  Shear  (N) 

(for  max.  t  =68.95MPa) 
xz 

2139.6 

1347.8 

653.9 

F±45o/Fs 

1.23 

1.31* 

1.35 

*The  60%  configuration  is  chosen  for  the  experiments. 
Note:  All  failure  stresses  are  for  uniaxial  tension. 

Lamina  Thickness  =  0.14  mm 
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Figure  1.  Minimechanics  Model  for  a  Notched  Laminate 


Figure  3. 


The  Various  Failure  Modes  in  a  Notched  Laminate 


Virgin  Data  Input 


transverse  Laminate  Failure  ? 


Figure  6. 


Flow  Chart  for  Failure  Mode 
and  Strength  Predictions. 


Distance 
from  Laminate 
Mid-plane 


0  =  45 


0  =  0 


aUn 

d  z 


0  =  -45 


Laminate 

Displacement’ 


Figure  7  .  Transverse  Displacement  Distributions  Near  Sp 
Edge,  [45/0/-45]  Laminate  (ref.  4  ) ,  Showing 

Average  Lamina  Displacements,  u  ,  First  Inter 
Shear  Strain,  y  .  and  Determination  of  Parame 


Elastic -Perfectly  Plastic  Approximation 


Figure  8. 


Different  Linear  Approximations  of  the  In-Plane  and 
Interlaminar  Shear  Strain  Curve 


Interlaminar 

plasticity 

between 

0°  and  ±45° 

layers 


Inplane  crack  in 
the  0°  layers  (C) 


Inplane  plasticity  in  the 
±45°  layers  ( ) 


Inplane  plasticity  in  the 
0°  and  90°  layers  (o^) 


Trial  Data  Input 


[0/145/90]^  Laminate 


a  (cm) 

Figure  10.  Variation  of  Applied  Laminate  Force  with 
the  Total  Axial  Damage  Zone  Size 


F 

(KN) 


Figure  12.  Variation  of  Applied  Laminate  Force  with 
the  Total  Axial  Damage  Zone  Size 


83 


(a  =■=  0.6  35  cm) _ 

Plasticity  in  the 
interface  between 
the  0°  and  ±45°  plies 

1 . . . . ^ 

Inplane 
plasticity 
in  the 

0°  plies 

tY  =  24.13  MPa 
xz 

[  0/0/ ±  4  5 ]  Laminate 
s 

_ ± _ 1 _ _  L  .  .  L__» 

Figure  13.  Variation  of  Applied  Laminate  Force 
with  Axial  Damage  Zone  Size 


Figure  17-  Fatigue  Failure  of  a  Notched  Boron/Epoxy 

shnwi na  Delamination  or  the  0°  Layers 


Laminate  Showing  Delamination  of 
from  the  Laminate  (S=.8)  (  ref.  2 


ftJMI] 


[0/+45/-45]  Laminate 
s 


0.281 


y  =  0.625 


Uniform  value  of  o  at  y  =  0.281 

yz 

from  the  approximate  analysis 


V  Uniform  value  of  at  y  =  0.625 

\/  from  the  approximate  analysis 


I  0\2 


-•—Edge 
of  the  notch 


0 .281 


-  0.625 


Figure  19  .  Average  a  for  Unit  Applied  Average  Stress 

y  z 

at  the  Interface  Between  0°  and  45°  Layers 
for  Different  Values  of  y. 


Element  Centroidal  ^yy  Values  in  the  Different 
Layers  for  Different  y  Values  from  the  Edge  of 
the  Hole 


Figure  2 0 . 


0.07 


0.06 


0 . 05 


0.04 


0.03 


0.02 


0.01 


Figure  2  3 . 


Circular  Hole 
x  =  0.125  (Edge  of  Hole) 
Order  of  Integration  =  4 


Variation  of  a  with  y  at  the  0/+45 
vz  J 

Interface  in  a  Notched  [0/+45/-45]  Laminate 


0.02 


Order  of  Integration  =  4 
x  =  0.125  (Notch  Edge) 


-0.04 


-0.06 


•0.08 


Slit  Notch 


-0.02 


Circular  Hole 


Figure  24.  Variation  of  o  with  y  at  the  +45/-45 

y  z 

Interface  in  Notched  [0/+45/-45]  Laminates 

b 

With  a  Circular  Hole  and  a  Slit  Notch 


Y 


Figure  26.  Variation  of  a  with  y  for  P  =  39.39E3 


c 


100 


101 


Figure  31.  "C"  Scans  for  a  [±45]  T-300/5208  Laminate  Panel 


Load  Cel  I 


Figure  32. 

Loading  Arrangement 
for  the  Modified 
IITRI  Compression 
Test  Fixture 


Actuator 
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106 


Figure  35.  [±45]  T-300/5208  Laminate  Failure  Modes  for  Uniaxial 

Tension  for  Two  Different  Specimen  Geometries 


Figure  36.  [±45^] g  T-300/5208  Laminate  Compression  Failures  Showing 

Microbuckling  of  Laminae 
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Figure  37.  [±45^g]  T-300/5208  Laminate  Failure  due  to  Surface 

Ply  Buckling 


Figure  38.  S-N  Data  for  [±252/9Q]  T-300/5208  Laminate 
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Figure  39.  [±259/90]  T-300/5208  Laminate  Edge  Failure  due  to 


tially  Notched  [±45/0~]  T-300/520 
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t 

Partial 
J  Notch 

+ 


Figure  41. 

C-Scans  of  Failed, 
Partially  Notched 

[ ±4  5/O3 ] s 
T-30  0/5208 
Laminates 


Interlaminar  Debond 


Implanted 
Crack 
In  0° 
Laminate 


Figure  42. 


Photomicrograph  of  a  Failed  [±45/0^] 
Partially  Notched  Laminate  (xlOO) 


113 


q  4f>  uej  4 g  |  Dnpisay 


116 


C  ycles 

Figure  45  .  Residual  Strength  -  N  Data  for  [0]  T-300/5208  Laminate 
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Figure  46.  [0]  T-300/5208  Laminate  Compression  Failures  during 

Residual  Strength  Tests 


T300/5208 

[0/90/Ok, 


o  §  X 
0  5>  O 

N 

o  " 

o  cr 


Figure  47  .  S-N  Data  for  [0/90/0] _  T-300/5208  Laminate 


iMDUdJiS  iDnpisay 


trength  -  N  Data  tor  [0/90/0] ~  T-300/5208  Laminate 


Figure  49.  [0/90/0] „  T-300/5208  Laminate  Failure  Modes  during 

£•  S 

Fatigue  (Inside  the  Tabs  for  S=0.8)  and  during 
Residual  Strength  Tests 
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